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FOREWORD 

This  report  to  the  Office  of  Naval  Research  Is  the  concluding  final 
report  for  a  two-year  theoretical  and  experimental  effort  carried  out  under 
Contract  N00014-80-C-0498,  and  covers  the  period  Hay  13,  1981  to  September  31, 
1982.  The  program  was  carried  out  at  the  Rockwell  International  Science 
Center  and  was  managed  by  Dr.  Paul  R.  Newman.  The  Principal  Investigators 
were  Hr.  Hark  D.  Ewbank  and  Dr.  Newman.  Prof.  Walter  A.  Harrison  of  Stanford 
University  was  a  consultant  on  the  theoretical  aspects  of  this  program.  Valu¬ 
able  contributions  were  also  made  by  Prof.  Eltan  Ehrenfreund  of  the  Technlon 
Institute  of  Israel,  and  Dr.  Pochl  Yeh,  Hr.  Randolph  L.  Hall  and 
Dr.  H.  Khoshnevlsan  of  the  Science  Center.  The  Contract  Honltor  for  the 
Office  of  Naval  Research  was  Dr.  George  Wright. 
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1.0  INTRODUCTION 


1.1  Program  Objectives 

The  long-range  goal  of  this  program  was  to  theoretically  Investigate 
the  microscopic  physical  origins  of  the  linear  dc  electro-optic  effect.  A 
secondary  objective  was  to  use  the  formalism  resulting  from  that  Investigation 
to  develop  a  general  methodology  with  predictive  capabilities  that  requires 
only  atomic  structure  and  elemental  composition  as  Input. 

Additionally,  as  new  Increasingly  high-performance  materials  are 
required  by  electro-optic  device  technology,  an  experimental  effort  was 
carried  out,  not  only  to  provide  feedback  for  the  theoretical  studies,  but  to 
provide  data  on  new  and  potentially  Interesting  Pockel's  electro-optic 
materials. 

1.2  Previous  Work 

The  previous  year's  work  was  first  devoted  to  separating  the  physical 
sources  of  the  electro-optic  effect  Into  those  which  result  primarily  from 
electronic  effects,  and  those  associated  with  lattice  dynamics  or  "ionlc- 
dlsplaclve"  effects.  The  Bond  Orbital  Model  *  was  then  modified  so  that  the 
electronic  contribution  to  the  electro-optic  effect  could  be  calculated  simply 
and  straightforwardly  using  "universal  atomic"  parameters  and  numerical  sums 
over  chemical  bonds  in  the  crystallographic  unit  cell.  This  theory  was  then 
tested  against  Te02  and  predicted  an  Identically  zero  electronic  contribution 
to  the  electro-optic  tensor  In  agreement  with  Klelnman's^  symmetry  relations. 

The  effort  next  focussed  on  the  ionic  dlsplaclve  part.  Here  a 
further  distinction  was  made  between  relative  motions  between  Ions  which  do 
not  result  In  a  change  In  the  unit  cell  dimensions  and  those  which  do.  The 
approach,  which  had  not  been  completely  evolved  at  the  conclusion  of  the  first 
year's  effort  was  to  seek  a  mechanism  to  calculate  the  microscopic  atomic 
"spring  constants"  and  relate  these  to  measurable  bulk  properties  such  as  the 
Reststrahl  frequency  or  elastic  constants. 
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Finally,  the  experimental  portion  of  the  program  was  directed  towards 
two  compounds:  Te02  and  Tl3AsSe3.  Measurement  on  Te02  yielded  the  original 
characterization^  of  the  linear  electro-optic  effect  In  terms  of  the  single 
non-zero  component  r41.  Preliminary  measurements  of  the  electro-optic  effect 
In  Tl3AsSe3  experienced  experimental  difficulties  associated  with  the 
semiconducting  nature  of  the  material. 

1.3  Accomplishments 

Significant  original  results  towards  understanding  the  electro-optic 
effect  have  been  made  during  this  period.  The  theoretical  work,  utilizing 
tight-binding  theory  and  Harrison's  "universal  atomic  parameters"1  has  suc¬ 
cessfully  yielded  a  formalism  for  quantitatively  calculating  both  the  elec¬ 
tronic  and  Ionic  contributions  to  the  electro-optic  second-order  suscepti¬ 
bility.  This  formalism  predicts  the  electro-optic  tensor  coefficients  util¬ 
izing  atomic  matrix  elements,  and  term  values,  together  with  crystallographic 
coordination  and  elemental  composition,  but  without  any  adjustable  parameters. 
As  was  Indicated  above,  the  electronic  contribution  is  calculated  in  a 
straightforward  manner  using  tight-binding  formalism  and  bond  sums  taken  over 
the  crystallographic  unit  cell.  Calculations  of  the  contributions  due  to  lat¬ 
tice  dynamics  has  also  been  successfully  reduced  to  practice  and  applied  to 
Te02.  The  formalism  involves  modelling  Interatomic  forces  (springs)  In  terms 
of  a  bond-stretching  force  and  a  bond-bending  force.  These  Interactomic 
forces  can  then  be  used  to  calculate  normal  acoustic  and  optical  lattice 
phonon  modes.  The  calculated  frequencies,  at  zone-center,  can  then  be  com¬ 
pared  with  data  derived  from  optical  measurements. 

Experimental  measurements  of  the  rc  electro-optic  tensor  component  of 
Tl3AsSe3  (TAS)  have  also  been  completed.  The  measurement  was  complicated  by 
the  semiconducting  nature  of  TAS.  When  several  kilovolts  were  applied  to  the 
sample,  appreciable  current  was  observed.  This  resulted  in  joule  heating  of 
the  sample  and  a  distortion  of  the  physical  dimensions  of  the  crystal  as  well 
as  concurrent  strain-optic  effects.  This  problem  was  eliminated  by  using  high 
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voltage  pulses  Instead  of  a  dc  voltage  and  yielded  the  "undamped"4  value  for 
rc- 

1.4  Technical  Issues 

Although  we  have  demonstrated  the  utility  of  our  approach  to  calcula¬ 
tions  of  the  electro-optic  effect  In  a  specific  case,  Te02,  the  generality  of 
the  method  remains  untested.  Further  work  on  a  variety  of  crystal  structures 
and  elemental  composition  needs  to  be  carried  out. 

1.5  Report  Summary 

The  remainder  of  this  report  Is  concerned  with  the  technical  details 
for  both  the  theoretical  and  experimental  aspects  of  this  program.  The 
material  contained  within  will  provide  the  basis  of  at  least  three  future 
publications,  which  are  currently  being  completed.  Finally,  a  brief  summary 
of  the  technical  progress,  along  with  some  projections  for  a  continuation  of 
work  in  this  general  area,  will  be  provided. 
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2.0  TECHNICAL  APPROACH 


2.1  Theoretical 

First,  the  connection  between  the  experimental  electro-optic  coeffi¬ 
cients  and  the  theoretical  second-order  susceptibility  tensor  is  developed. 
Then  a  detailed  formulation  for  the  dielectric  susceptibilities  Is  made  in 
terms  of  tight-binding  theory.  Initially,  the  linear  susceptibility  is 
treated  and  subsequently  the  model  is  extended  to  higher-order  suscepti¬ 
bilities.  In  particular,  the  electro-optic  susceptibility  is  examined  and  a 
technique  for  evaluating  the  lattice  contribution  is  presented.  Finally,  a 
specific  application  of  the  theory  is  carried  out  on  tellurium  dioxide. 


2.1.1  Second-Order  Susceptibilities  and  Electro-Optic  Coefficients 

Experimentally,  the  determination  of  the  change  in  refractive  index 
of  a  crystal  due  to  an  applied  dc  electric  field  provides  a  mechanism  of 
measurement  of  the  electro-optic  effect.  The  electro-optic  coefficient,  r^jk, 
Is  defined  by  the  relation4-*! 

A(B),J  -  I  rtjkEk  (1) 

where  AB^.  is  the  change  in  the  relative  optical  dielectric  impermeability  and 
Ek  is  the  applied  electric  field.  The  impermeability  Is  related  to  the 
Inverse  dielectric  tensor,  e"*,  and  the  refractive  index,  n,  by  the  expression 


Bu  ■  <*'%  ■ 


.-2> 


(2) 


A  change  in  impermeability  can  be  expressed  in  terms  of  an  index  change  as 


AB  =  -2An/n3 


(3) 
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and  explicit  relationships  between  in  and  ABjj  have  been  derived  for  optically 
isotropic,  uniaxial  and  biaxial  crystals  with  the  light  propagating  In  an 
arbitrary  direction.6 

The  principal  dielectric  constant  or  refractive  index  can  In  turn  be 
related  to  the  susceptibility,  x»  as 

n  =  /e  =  /I  +  4*x  (4) 

in  CGS  units.  To  define  higher-order  susceptibilities,  the  induced 
polarization,  P^,  is  expanded  in  terms  of  increasing  powers  of  electric  field1 
as 


Pi  =  l  xij>Ej  +  ^  xijk  EjEk  +  xi  jkl  Ej  EkEl  + 


,(?> 


J3) 


(5) 


where  xjj^  Is  the  linear  or  first-order  susceptibility,  xjjk  is  the  second- 
order  susceptibility,  etc.  This  form  does,  however,  neglect  any  frequency 
dependence.  More  specifically,  the  induced,  second-order  polarization  at  any 
given  frequency  Is  a  function  of  all  three  wave-mixing  combinations  which 
contribute  to  that  frequency,  i.e..  Integrate  the  product  of  the  second-order 
susceptibility  with  two  electric  fields  over  frequency.7  By  considering  three 
special  cases  in  the  frequency  domain,  three  distinct  second-order  suscepti¬ 
bilities  have  been  defined:6 


=  l  x®  (2<o)Ej(<»)Ek(ta>) 
J  * 


(6) 


0)  =  l  xjjk  (0)Ej{ui)Ek(u)) 
J  k 


(7) 
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=  l  xjkju)  Ej ( so)  Ek(0)  (8) 

J 

which  correspond  to  second  harmonic  generation,  optical  rectification  and  the 

linear  electro-optic  effect,  respectively.  The  electro-optic,  second-order 

susceptibility,  as  defined  in  Eq.  (8),  must  be  symmetric  in  he  indices  i  and  j 

in  order  that  the  total  optical  susceptibility  tensor,  +  1  x-l-  El.  remain 

^  J  £  1 X  J  K 

symmetric.  This  contrasts  with  the  second  harmonic  generation  and  optical 
rectification  susceptibilities  which  are  symmetric  in  all  three  indices  (i ,  j 
and  k).  Keeping  these  symmetry  relations  straight  is  important  when  dealing 
with  "condensed"  notation^  for  the  indices  because  the  condensing  always 
occurs  on  the  two  symmetric  indices.  We  shall  be  concentrating  on  the 
electro-optic  second-order  susceptibility,  which  can  be  directly  related  to 
the  electro-optic  coefficient, 7«9*10  in  CGS  units,  by 


rijk  =  4’4k]/en/ejj 


(9) 


1  +  is  an  optical  dielectric  constant.  Note  that  both 

are  symmei 

dielectric  constants  appearing  in  Eq.  (9).  The  theory  developed  below  derives 
expressions 
experiment. 


2 

where  e . .  =  n . 

M  /£  \  ' 

r-jjk  and  x^jl10)  are  symmetric  in  indices  i  and  j  which  correspond  to  the 


expressions  for  xjj^(“)  and  therefore  Eq.  (9)  will  serve  as  the  connection  to 


2.1.2  Tight-Binding  Calculations  of  the  First-Order  Susceptibility 

For  nominal  optical  electric  fields,  the  linear  susceptibility  is 
usually  related  to  the  interband  absorptions.  These  ehsorptions  contribute  to 
an  imaginary  component  of  the  susceptibility  which  corresponds  to  a  polariza¬ 
tion  density  that  is  out  of  phase  with  the  driving  electric  field.  This 
absorptive  component  can  be  written  as 
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Im 


* 

L*ij  J 


IT 

15 


s  ■Hrl*v(k*»<*ll(if)l(|3)  ^-l»c(M> 


v,c,kk' 


x  6(Ec(i< 1 )  -  Ey(it)  -  hu>) 


(10) 


where  e  is  the  electronic  charge,  m  is  the  mass  of  electron  and  h  is  Planck's 
constant  divided  by  2 it.  The  dipole  matrix  element  is  represented  by  the 
gradient  operator,  $,  coupling  the  occupied  (valence)  wave  function, 

|*v(ic)>,  with  the  unoccupied  (conduction)  wave  function,  (£')>.  The  delta- 
function  gives  a  contribution  only  when  the  photon  energy  (hu>)  equals  the 
energy  difference  between  the  coupled  states.  Finally,  the  summation  Is  over 
wave-vectors  of  all  valence  and  conduction  states  in  a  volume,  n. 

The  real  part  of  the  linear  susceptibility  is  related  to  the 
imaginary  component  via  the  Kramers-Kronig  relations.1  By  assuming  that  the 
dipole  matrix  elements  are  independent  of  frequency,  the  real  part  of  *(!)  can 
be  written  as 


(1)  s  4eV 

*1  j  2 

J  nrn 


<*(*)  l-drl \-$q\ \(*)> 

,.c.f  <Ec<k*>  -  E,ff)  )J 


(11) 


in  the  limit  that  the  photon  frequency  is  much  less  than  the  energy  difference 
between  coupled  states  or  zero-frequency  limit,  (i.e.,  hu>  is  much  less  than 
the  lowest  absorption  energy  but  still  greater  than  the  Reststrahl  vibrational 
energies).  An  extra  factor  of  two  has  been  Included  In  Eq.  (11)  to  account 
for  the  two  spins  associated  with  each  state  and  the  summation  over  wave- 
vector,  ft',  has  been  omitted  by  considering  only  "vertical  transitions." 

Also,  note  that  the  designation  of  the  real  part  is  not  explicitly  shown  in 
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Eq.  (11)  since  only  the  real  components  of  the  susceptibilities  will  be 
considered  in  the  rest  of  the  formalism  that  follows. 

2. 1.2.1  Tight-Binding  Theory 

In  the  context  of  tight-binding  theory,  the  dipole  matrix  elements  In 
Eq.  (11)  can  be  expanded  In  terms  of  localized  Wannier  functions,  which  in 
turn  are  written  as  a  summation  of  Bloch  functions.  Specifically,  the 
valence-band  state,  \py ( it) ,  can  be  expressed  as 


11c. It 

Uv(k+)>  =  /jf  m  ,av(^)>e  *m  (12) 

where  av((fm)  is  the  valence  Wannier  function  located  at  site,  in  the  m-th 
primitive  cell  and  the  summation  is  over  N  primitive  cells.  The  conduction- 
band  wave  function  is  represented  in  an  analogous  fashion.  By  carrying  out 
this  expansion  of  the  dipole  matrix  elements  in  Wannier  functions,  the  double 
summation  over  sites  in  each  of  N  primitive  cells  reduces  to  N  times  a  single 
summation.  This  is  because  for  any  particular  dipole  matrix  element  involving 
two  Wannier  functions  there  are  only  N  identical  matrix  elements  of  Wannier 
functions  on  nearest  neighbor  sites  which  are  non-negligible.  Then,  the 
summation  over  N  primitive  cells  becomes  a  summation  over  nearest  neighbors. 
Finally,  the  dipole  matrix  element  in  Eq.  (11)  is  written  in  terms  of  Wannier 
functions  as 


.  a  +  (nn)  *  a  * 

<*c<k>liri*v<k)>  •  I  V*i>lisi!rlaA»e 

1m  1 


(13) 


where  ^  specifies  a  "central"  primitive  cell  conduction-state  site  and  the 
summation  is  over  only  those  valence-state  sites  ^  which  are  nearest  neigh¬ 
bors  (nn)  to  this  conduction-state  site  ^  (including  the  valence  states  on 
the  same  site). 
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Furthermore,  when  this  dipole  matrix  element  Is  substituted  back  into 
the  expression  for  the  linear  susceptibility,  it  Is  recognized  that  the  wave- 
vector  dependence  of  the  energy  difference  denominator  is  negligible  (since 
the  energy  variation  with  if  Is  small  compared  to  the  energy  separation  between 
valence  and  conduction  states).  Also,  the  summation  over  wave-vector  can  be 
eliminated  because  the  sinusoidal  factor  Integrated  over  if  contributes  only  at 
a  single,  localized  valence-state  site.  Then  the  linear  susceptibility  can  be 
expressed  as 


where  the  summation  includes  only  those  nearest  neighbor  valence  and 
conduction  Wannier  functions,  located  at  sites  and  ,  respectively. 

Typically,  the  Wannier  functions  have  been  chosen  by  incorporating 
nature  of  the  chemical  bonding  of  the  particular  crystal  under 
consideration.  For  example,  the  tetrahedral  semiconductors  have  four  bond 
orbitals  per  primitive  cell  from  which  four  valence-band  Wannier  functions  are 
derived  from  the  four  antibond  orbitals.  By  focusing  on  only  those  dipole 
matrix  elements  of  the  gradient  between  bonding  orbitals,  |b>,  and  anti  bonding 
orbitals,  |a>,  in  the  same  bond,  the  summations  over  nearest  neighbor  valence 
and  conduction  states  became  a  summation  over  bonds.  For  ionic  compounds, 
such  as  the  alkali  halides,  the  valence-  and  conduction-band  Wannier  functions 
are  adequately  represented  by  the  hallde-p  and  alkall-s  atomic  orbitals, 
respectively.  In  this  case,  the  summation  over  band  states  became  a  summation 
of  spo  couplings  over  the  six  nearest  neighbors. 

However,  it  is  not  necessary  to  specify  the  chemical  bonding  of  the 
crystal  at  this  point  in  the  formalism.  Nothing  Is  lost  In  generality  by 
representing  the  Wannier  functions  by  a  complete  set  of  valence-electron 
atomic  orbitals:  |s>,  |px>,  |py>,  |pz>,  etc.  For  example,  in  the  above  case 
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of  the  tetrahedral  semiconductors,  hybrid  SP3  bonding  and  antibonding  orbitals 
(and  corresponding  matrix  elements)  are  simply  derived  from  the  appropriate 
linear  combinations  of  atomic  orbitals.  Then,  from  Eq.  (14),  the  dipole 
matrix  elements  are  between  atomic  orbitals  and  the  summation  over  nearest 
neighbor  valence  and  conduction  Wannler  functions  becomes  a  sum  over  all 
nearest  neighbor  pairs  of  atomic  orbitals  in  the  primitive  cell. 

2. 1.2. 2  Simplified  Evaluation  of  the  Dipole  Matrix  Element 

The  dipole  matrix  element  involving  atomic  orbitals,  |a>  and  |b>,  on 
neighboring  atoms  can  be  approximated  In  a  very  simple  form  in  terms  of 
known1**1  parameters.  These  known  parameters  are  the  atomic  term  values  (ea 
and  eg,  corresponding  to  each  atomic  orbital)  and  the  matrix  element  of  the 
crystal  Hamiltonian  between  the  two  atomic  orbitals  (Vag  *  <a|H|B>).  In 
addition,  the  two  atomic  orbitals  are  assumed  to  be  orthogonal:  <a|&>  =  0. 

To  obtain  the  desired  dipole  matrix  element,  the  eigenstates  in  the  LCAO 
context  of  this  isolated  pair  of  atomic  orbitals  will  first  be  constructed. 
Then,  the  molecular  dipole  matrix  will  be  related  back  to  the  dipole  matrix 
element  between  two  atomic  orbitals. 

The  eigenvalues  are  simply  given  by 


=  e6)/2  *  /lEa  '  V 


where  the  subscript  c  and  v  are  refer  to  the  conduction  and  valence,  respec¬ 
tively.  The  eigenstates  are  written 


|v>  *  lT|o>  +  U'|6> 


|c>  -  UCa(o>  +  Ugl  8> 
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and  the  LCAO  coefficients  are  given  by 


Uva .  +UC6  -  /rrr-y72 

uv0  -  -uca  -  /TT-r-^W  (17) 


where 


%  •  -  s>2  ♦  ']'1/2  »8> 

and  represents  the  "polarity"  of  the  coupling. 

Making  use  of  a  general  relation12  among  matrix  elements  between 
eigenstates. 


<vbfHC>  =  -  ey)<v|xi|c>  (19) 

the  dipole  matrix  element  with  a  gradient  operator  Is  replaced  by  a  matrix 
element  with  a  position  operator.  The  position  matrix  element  of  the  two 
eigenstates  is  immediately  expressed  as 


<v|x1|c>  =  /{ 1  +  Op){l  -  ap)  (<g|x1|6>  -  <a|xi|a>)/2  +  <^)<o|xi|0>  (20) 


where  |a>  and  | 0>  are  assumed  to  be  real  which  implies  that  <e|x^|a>  = 

<a|x^ | 6>.  This  last  term  vanishes  if  the  atomic  energies  of  the  states  are 
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the  same  (l.e.,  ap  =  0)  or  If  the  two  states  are  symmetric  about  the  midpoint 
between  them.  (Note  that  since  the  two  atomic  states  are  orthogonal,  this 
matrix  element  in  the  last  term  is  independent  of  the  origin  of  coordi¬ 
nates.)  As  in  early  treatments,*  this  term  is  assumed  to  be  zero  (the  only 
step  In  the  derivation  of  the  dipole  matrix  element  of  two  atomic  orbitals 
that  Is  not  rigorously  justified).  Then  the  first  term  In  Eq.  (20)  Is 
proportional  to  the  i-th  component  of  the  interatomic  distance  vector,  <f, 

d.  =  <sjx1|S>  -  <o|xiia>  (21) 

permitting  Eq.  (20)  to  be  written  as 

<v|xi|c>  =  d./2  (22) 

The  gradient  matrix  element  of  two  eigenstates  in  Eq.  (19)  is  expanded,  in  a 
similar  fashion  to  Eq.  (20),  as 

<V|_L|0  .  /rr^  «e|-S-|e>  -  <«l^-|o»/2  ♦  (23) 

where  a  partial  integration  has  been  used  to  show  that  <e|3/3x^|a>  ■ 
-<a|3/3xi|e>  In  the  last  term.  Additionally,  the  same  partial  Integration 
shows  that  both  <a|3/3x^|a>  and  <s|3/3x^|8>  vanish,  so  Eq.  (23)  becomes 

.  <0|_|-|6>  .  (2„ 
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Substituting  Eqs.  (22)  and  (23)  back  Into  Eq.  (19)  and  solving  for  the  dipole 
matrix  element  between  atomic  states  gives 


<a|~-|B>  •  (m/«2)(ec  -  ev)/m?d1/2 


(25) 


This  expression  further  simplifies  by  using  Eqs.  (15)  and  (18)  to  give 
(ec  -  E^,)  =  2VagA‘  -  Op  which  yields 


(26) 


for  the  final  form  of  the  dipole  matrix  element  involving  two  atomic  orbitals. 


2. 1.2.3  Simplified  Linear  Susceptibility  Formula 

This  form  for  the  dipole  matrix  elements  (Eq.  (26))  can  be  utilized 
In  the  tight-binding  formulation  of  the  first-order  susceptibility  (Eq.  (14)) 
when  atomic  orbitals  are  substituted  for  the  Wannier  functions.  The  energies 
in  the  denominator  of  Eq.  (14)  can  be  adequately  represented  by  the  eigen¬ 
values  of  a  pair  of  coupled  states  In  Eq.  (15)  and  the  resulting  linear  sus¬ 
ceptibility  is  written 


(1).  «ji  <"r">  v.e  Vi 

«•»  c(.8  -  gz  * 


where  the  summation  is  over  nearest  neighbor  atomic  orbitals. 


(27) 
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2. 1.2. 4  Applications 

This  formula  Is  directly  applicable  to  any  specific  material  with  all 
the  necessary  parameters  known  once  the  crystal  geometry  has  been  determined. 
The  only  complication  Is  the  possibility  of  a  large  number  of  terms  entering 
the  summation  over  atomic  orbitals.  Even  further  simplifications  can  result 
If  assumptions  about  the  electronic  band  structure  for  the  particular 
compounds  of  interest  are  incorporated  at  this  point. 

Consider,  for  example,  the  alkali-halides.  First,  note  that  the 
covalent  energy  Is  always  much  less  than  the  polar  energy  for  all  pairs  of  s 
and  p  states.  Consequently,  the  Va@  term  in  the  denominator  of  Eq.  (27)  Is 
negligible.  Also,  the  only  Important  atomic  states  which  enter  are  the  s- 
state  on  the  alkali  atom  and  the  p-state  on  the  halide  atom.  This  is  because 
the  energy  difference  between  other  pairs  of  nearest  neighbor  atomic  states  Is 
significantly  larger.  Finally,  in  a  principle  axis  coordinate  system,  the 
Interatomic  distance  vectors  are  directed  along  one  of  the  principle  axes 
which  means  that  d^dj  =  d2  -  a2/4  only  when  1  -  j.  Since  there  are  six 
nearest  neighbors  surrounding  each  of  the  two  atoms  in  the  primitive  cell. 


.2.4  2 
4e  *  nsp o 


m2d5(es  - 


(28) 


where  n  =  2d3  has  been  used  and  numerical  values  for  the  coefficient,  n^, 
are  given  In  Ref.  1.  This  prediction  has  been  compared  to  experiment 
elsewhere.* 

The  case  of  the  binary  semiconductors  Is  almost  as  simple.  Here, 
hybrid  orbitals  are  formed  and  the  Interatomic  matrix  element,  which  Is 
replaced  by  the  covalent  energy  (V2),  is  no  longer  negligible  and  gives 


/Te2V2 

64d(V2  +  V2)3'2 
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with  fl  -  2a^  =  128  d3/3/T  for  four  hybrid  bonds  and  where  V2  *  -3.22  fi^/fmd^) 
and  V3  =  (e^  +  3ep  -  e|  -  3ep)/8  are  the  covalent  and  polar  energies, 
respectively,  as  previously  defined. ^ 

2. 1.2. 5  Dipole  Moments 

Harrison  has  utilized  the  concept  of  a  dipole  moment  of  a  chemical 
bond  when  considering  the  susceptibility  of  tetrahedral  semiconductors.  The 
1th  vector  component  of  this  dipole  moment  has  been  written  as 

P1  a  ed1ap  (30) 

where  ap  =  +  V|  is  the  polarity  of  the  bond.  Then,  the  total  induced 

polarization  density  Is  just  the  summation  over  bonds  of  each  dipole  moment 
averaged  over  the  volume  occupied  by  the  bonds.  This  dipole  moment  concept 
can  also  be  employed  within  the  context  of  the  above  formalism.  Here,  the 
"dipole  moment"  represents  the  Induced  coupling  between  two  atomic  orbitals 
and  can  be  written  as  In  Eq.  (30)  but  with  the  "polarity"  given  by  Eq.  (18). 
The  connection  between  the  bond  dipole  moment  and  the  susceptibility  can  be 
established  rigorously  as  follows.  As  before,  the  total  Induced  polarization 
density  is 


P1  "  Z  J  p1 


where  fl  Is  the  volume  occupied  by  the  coupled  states.  Now,  with  reference  to 
Eq.  (5),  each  dipole  moment  Is  written  as  a  Taylor’s  series  expansion  of  the 
electric  field  and  Eq.  (31)  becomes 


2 

1  1  ^  P. 
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where  3pj/3Ej  Is  the  1-th  component  of  the  contribution  to  the  first-order 
susceptibility  from  each  dipole  moment.  After  Interchanging  the  summation 
over  coupled  states  with  summations  over  vector  components  In  Eq.  (32)  and 
comparing  back  to  Eq.  (5),  the  first-order  susceptibility  Is  then  Identified 
as 


m  .ij.ii 

xu  a 


cs*T 


(33) 


To  verify  this  expression,  one  can  carry  out  the  dlfferentltatlon  according  to 
the  chain-rule: 


where  the  Independent  variables  are  d^ ,  Va5  ~  1/d2  and  eoe  =  (ea-  eg)*  If 
the  frequency  of  the  an  optical  electric  field  Is  above  the  Reststrahl 
vibrational  frequencies,  the  atoms  do  not  respond  and  hence  the  change  In  the 
m-th  component  of  the  Interactomlc  distance  vector  with  respect  to  this 
optical  field,  (3dm/3Ej),  Is  zero.  Consequently,  only  the  last  term  In  Eq. 
(34)  Is  nonzero.  Note  that  the  energy  difference,  ea&.  Is  modified  by  the 
dipole  energy,  -ecf*!f,  when  the  optical  field  Is  applied.  This  last  term  Is 
easily  evaluated  and  the  resulting  linear  susceptibility  Is,  of  course, 
Identical  to  Eq.  (27). 
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2.1.3  Higher-Order  Susceptibilities 

The  various  second-  and  third-order  susceptibilities  can  be  treated 
In  a  similar  fashion  to  the  linear  susceptibility  in  Eq.  (33).  As  a  note  of 
caution,  the  particular  Taylor's  series  expansion  (as  used  in  Eq.  (32))  must 
reflect  the  experiment  being  considered.  For  example,  in  second  harmonic 
generation,  there  is  only  one  electric  field  which  Is  the  optical  field  and 
the  Taylor's  series  expansion  in  three  independent  variables  (the  three 
components  of  this  optical  field)  as  written  In  Eq.  (32)  Is  appropriate. 
Then,  the  second  harmonic  susceptibility  is  immediately  recognized  as 


SHG  _  1 
xijk  ”  2fl 


l 

cs 


a2P. 


3EjSEk 


(35) 


where  the  full  tensor  notation  has  been  explicitly  retained.  Again,  perform¬ 
ing  the  differentiation  Is  relatively  trivial  when  the  change  in  Interatomic 
distance  vector  is  zero,  as  it  is  for  the  case  of  second  harmonic  genera¬ 
tion.  The  result  can  be  written 


SHG  6e^ 

xijk  -  n  J 


(eB' 


e  )\r  .d^d. 
a'  o0  1  J  k 

~T - 


CK 


t  ) 
a 


(36) 


and  Is  obviously  symmetric  in  Indices,  i,  j  and  k  which  is  consistent  with 
Kleinman's  symmetry  condition.2 

2.1.4  Electro-Optic  Susceptibility 

As  another  example,  consider  the  linear  dc  (Pockel's)  effect.  The 
electro-optic  susceptibility  Is  related  to  the  changes  Induced  in  each  dipole 
moment  by  both  the  optical  and  dc  electric  fields.  This  can  be  seen  mathema¬ 
tically  by  writing  out  the  Taylor's  series  expansion  for  the  i-th  vector 
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component  of  the  dipole  moment  In  terms  of  the  six  Independent  variables 
(which  correspond  to  the  three  components  of  both  the  optical  and  dc  electric 
fields): 


p1<r,pt.fdc> 


p,  (f°pt  -  o.r*  -  o)  !  )Ejopt  *  i 

1  M  »E°pt  J  j=l 


(■ 


3pi 


J 


+ 


r.2 

3  p. 


‘5E?pt»E 


opt 

k 


+ 


3  3 

Ml)  I  I 

j=l  k=l 


,2n 

3  Pi 


3E0pt3E°Pt 

J  K 


)E°ptE^c  +  ... 


(37) 


The  contributions  from  each  dipole  moment  to  the  electro-optic  susceptibility 
arise  from  terms  proportional  to  E°ptE^C,  i.e.,  the  last  summation  in  Eq.  (37). 
Consequently,  the  electro-optic  susceptibility  can  be  written 


E0  1  r  _ _ 

Xikj  ‘  n  cs  3E?pt 


(38) 


where  the  summation  is  over  all  nearest-neighbor  coupled  atomic  states.  The 
derivative  in  Eq.  (38)  can  be  expressed  in  terms  of  a  set  of  known  parameters 
by  carrying  out  the  differentiation  on  Eq.  (30)  as  follows.  The  1-th  com¬ 
ponent  of  the  dipole  moment  in  Eq.  (30)  is  a  function  of  the  variables  di,  Va0 
and  tao.  We  apply  the  chain  rule  to  convert  the  derivative  with  respect 
to  E°pt  into  separate  derivatives  with  respect  to  these  three  variables. 

Noting  that  d^  and  VaB  are  constant  in  the  presence  of  a  high  frequency 
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optical  field,  E?^,  we  find  that  the  optically-induced  change  In  the  dipole 
moment  is1 


3Pj 

3E9P1 

J 


■  ( 


as?'1* 


)  - 


(e2  + 

*  aB 


<s>3/2 


(39) 


This  represents  the  contribution  of  each  dipole  moment  to  the  first-order 
susceptibility.  By  contrast,  the  dc  electric  field  induces  an  additional 
change  in  the  dipole  moment.  Involving  terms  in  which  include  the  lattice 
response.  Using  the  chain  rule  with  four  variables  (di ,  dj,  Vag  and  ea0), 
this  second-order  effect  becomes 


3di  3 

( — jr)  '3d- 
3Edc  3d. 


3di  a  3V„« 

+  r _ L-)  __L  +  r _ 2§i 

dcJ  3d.  1 


3E 


3E 


dc  ;  3V 


a0 


3c 


+  i-g) 


3Ek  3E“3 


4ezd1dj 

J 

(ea8  +  4VaB^ 


-2V  .  3  3d 

_ C*B  V  A  !  m 

^2  f.  dm  ^  Fdc 
d  m=l  3E  ^ 


(40) 


The  evaluation  of  the  lattice  response  to  the  dc  electric  field  (i.e.,  the 
factors,  3dn/3Edc  will  be  treated  as  a  separate  problem  in  Section  2.1.5. 

Then,  the  remaining  derivatives  in  Eq.  (40)  can  be  expressed  entirely  in  terms 
of  known  parameters,  assuming  the  crystal  structure  has  been  specified. 
Performing  these  differentiations  yields 


cn  aJ  (nn)  V2  (  3  3d  r  , 

x1kj  =  —  l  — i - — 57?  ! 3edi didkEa6  +  ^  ( — 3cH^e' 

1kJ  n  a,8  (c2e  +  4V26)5/Z  (  J  aB  m=l  3eJC  L 


2  +  4V2  ) 

a0  aB7 


<v v ♦  Vi»>  *  4<2vl8  -  '>idjVdi! 


(41) 
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2 

where  5^  is  the  Kronecker  delta,  eag  =  (ea  -  eg)  and  the  first  summation  is 
over  coupled  nearest-neighbor  atomic  orbitals. 

2.1.5  Lattice  Contribution  to  the  Electro-Optic  Effect 

First,  the  dynamical  matrix  method  will  be  used  to  describe  the  lat¬ 
tice  equations  of  motion  in  three  dimensions  for  a  generalized  crystal  struc¬ 
ture.  This  method  will  employ  only  nearest-neighbor  bond-stretching  and  bond- 
bending  force  constants*  for  simplicity.  Only  the  phonon  modes  at  zero  wave- 
vector  will  enter  the  calculation  of  the  electro-optic  effect;  however  the 
entire  lattice  spectrum  is  investigated  in  order  to  test  this  force-constant 
model.  Then  the  force  on  each  atom  in  the  unit  cell,  due  to  the  dc  electric 
field,  will  be  quantified  by  assigning  an  effective  charge  to  each  atom. 

Hence,  the  displacement,  in  response  to  the  dc  electric  field,  of  every  atom 
in  the  unit  cell  will  be  determined. 

2. 1.5.1  Lattice  Equations  of  Motion 

Consider  a  crystal  consisting  of  N  atoms  per  unit  cell.  When  the 
lattice  is  perturbed  (by  a  dc  electric  field,  for  example),  each  atom  in  the 
unit  cell  will  undergo  elastic  motion  in  response  to  this  perturbation.  This 
elastic  motion  can  be  modelled  using  a  nearest-neighbor  formalism  similar  to 
that  used  in  Eq.  (41).  In  this  model  each  atom  in  the  unit  cell  experience 
two  types  of  forces  when  it  is  displaced  from  its  equilibrium  position.  The 
first  is  a  bond-stretching  force  and  is  simply  due  to  the  radial  change  in 
interatomic  distance  between  the  m-th  atom  in  the  unit  cell  and  each  of  its 
nearest  neighbors  (atom  l),  as  shown  in  Fig.  la.  The  second  type  force  is  a 
bond-bending  force  and  is  due  to  the  change  in  bond  angles.  For  each  atom  in 
the  unit  cell,  there  are  two  different  kinds  of  bond  angles  to  consider;  one 
Involves  two  first-neighbors  (atoms  i  and  n)  while  the  other  Involves  a  first- 
nearest  neighbor  (atom  n)  and  a  second-nearest  neighbor  (atom  m').  These  two 
contributions  to  the  elastic  motion  of  the  m-th  atom  in  the  unit  cell  are 
Illustrated  in  Figs,  lb  and  lc,  respectively. 
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Equilibrium  and  non-equilibrium  atomic  configurations  of  nearest 
neighbor  atoms  for:  (a)  bond- stretching  forces,  (b)  bond-bending 
forces  with  two  first-nearest  neighbors,  and  (c)  bond-bending  forces 
with  a  first-  and  a  second-nearest  neighbor. 
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The  potential  energy  associated  with  the  bond-stretching  between 
atoms  n  and  i  can  be  written 

■  7  co«  <42> 

where  Cq  is  defined1  as  the  radial  force  constant,  Ad^  is  the  change  in 
magnitude  of  the  interatomic  distance  and  d^  is  the  equilibrium  interatomic 
distance.  Similarly,  the  potential  energy  associated  with  the  bond-bending 
between  atom  m  and  nearest  neighbor  atoms  i  and  n  is 

6«nn  ■  7  V9L  (43> 

where  Cj  is  defined1  as  the  angular  force  constant  and,  ab^  is  the  change  in 
the  bond  angle.  Finally,  the  potential  energy  associated  with  the  bond¬ 
bending  between  atom  m,  first-nearest  neighbor  atom  n  and  second-nearest 
neighbor  atom  m'  can  be  written  in  a  fashion  analogous  to  Eq.  (43). 

Using  these  three  contributions  to  the  potential  energy,  the  total 
force  on  the  m-th  atom  in  the  unit  cell  becomes 


L  ■  - I  -i  :  '  - 1  i  *■ 


l  n *l 


n  m'*l 


7  l  -  7  LL  (Cl>»n?<i9lmn> 


l  T\*i 


7  J  .X  <Cl>mnm-5!49L') 


n  m  m 
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where  the  summation  for  the  bond-stretching  contribution  is  over  all  nearest 

neighbor  atoms  (t),  the  summation  for  the  first  bond-bending  contribution  is 

over  all  combinations  of  nearest-neighbor  pairs  (*  and  n)  and  the  summation 

for  the  second  bond-bending  contribution  includes  each  first-nearest  neighbor 

atom  (n)  paired  with  all  of  its  first-nearest-neighbor  atoms  (m')  excluding 

atom  m.  The  gradient  in  Eq.  (44)  is  with  respect  to  the  displacement  vector 

associated  with  atom  m.  Note  that  there  is  only  one  set  of  angular  force 

constants  for  (Cj)^  and  (Ci)mnm^.  The  differentiations  indicated  in 

Eq.  (44)  can  be  carried  out  in  a  straight-forward  manner.  First,  for  the 

bond-stretching  contribution,  the  change  in  magnitude  of  the  interatomic 

distance  can  be  written  explicitly  as  a  function  of  the  displacement  vector, 

u  ,  of  atom  m: 
m 


^ 9m  ^«m  ’  -  ’  u2m^d«m 


where  a  binomial  expansion  has  been  used  and  only  terms  up  to  first  order  in 
the  displacements  were  kept.  Here  and  throughout  a  short-hand  notation  is 
used  for  the  difference  between  two  vectors,  e.g.,  the  displacement  vectors  in 


Eq.  (45)  enter  as  Uj. 
expressed  as 


2  “*  -  v 


Then  the  gradient  of  Ad...  is  simply 


~  "2^*m^*m  '  u*m^d*m 


The  bond-bending  contributions  are  considerably  more  complicated,  mainly  due 
to  the  algebra.  The  change  In  bond  angle  can  be  related  to  the  displacement 
vectors,  um,  as  follows.  In  the  equilibrium  configuration,  the  cosine  of  the 
bonding  angle  is 
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cos  "wn  = 


written  in  terms  of  the  interatomic  distance  vectors.  For  the  perturbed 
configuration,  the  cosine  of  the  bonding  angle  becomes 


cos  9'  =  J  (l  -  •  u„  /q  -  cf  •  u  la  ) 

2mn  L  in  nm'  2m  2m'  2m  nm  nm'  nm' 


+  S  •  u  +  J  •  u  ]/(d  _d  ) 
2m  nm  nm  2m J/  v  2m  nm' 


again  using  the  binomial  expansion  and  discarding  terms  second  order  or  higher 
in  the  displacements.  The  small  change  in  bond  angle  is  related  to  the 
difference  in  the  two  cosines  (Eqs.  (47)  and  (48))  via  the  definition  of  the 
derivative  for  the  cosine.  This  change  in  angle  can  then  be  written  as 


A02mn  “  ‘^nm  "  ^2m  cos  02mn^  "  u2n/c*2m  +  ^2m  ”  ^nm  cos  02mn^ 


*  u^d  1/sin  9. 
nm  nm  2mn 


where  d„  and  d  are  unit  vectors  directed  between  nearest-neighbor  atoms. 

2m  nm 

Finally,  the  two  gradients  with  respect  to  the  displacement  vector,  um,  of  the 
square  of  the  change  in  bond  angle  (after  a  fair  amount  of  algebra)  becomes 
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(50) 


for  the  case  of  two  first-nearest  neighbors  to  atom  m  and 


^A9mnm^  =  "2J  '  2 


Sl" 

mnm  mn  m  n 


d2  d 

.n.  ^  COS2  0  T~~  &  '  cos  9 

^  mn  mnm  d  mn  mnm 
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(51) 
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for  the  case  of  one  first-nearest  neighbor  and  one  second-nearest  neighbor  to 
atom  m. 

By  substituting  Eqs.  (46),  (50),  and  (51)  back  into  Eq.  (44),  the 
generalized  equation  of  motion  for  the  m-th  atom  in  the  unit  cell  is  obtained. 
While  this  equation  of  motion  appears  to  be  rather  complicated  on  the  surface, 
closer  inspection  reveals  that  most  of  the  factors  (e.g.,  and  djp,,  9^  and 
0mnm')  are  direct1y  related  to  the  unperturbed  crystal  geometry  and  also 
that  this  equation  is  nonlinear  in  the  displacements.  The  only  parameters 
which  are  not  immediately  in  hand  are  the  set  of  bond-stretching  and  bond- 
bending  force  constants  ((Cq)^  and  (C^)^  respectively).  For  simple  crys¬ 
tals,  such  as  the  binary  semiconductors,  the  bonds  are  all  identical  (or 
nearly  so)  and  only  a  single  pair  of  numbers  is  needed  for  C0  and  Cj.  In  more 
complex  crystal  structures,  the  bonds  are  no  longer  equivalent  and  the 
degeneracy  in  the  force  constants  is  broken.  Methods  for  determining  these 
force  constants  have  been,  and  are  still  being,  investigated.13  Consequently, 
for  the  problem  being  addressed  here,  it  is  assumed  that  these  force  constants 
are  known. 


2. 1.5. 2  Lattice  Dispersion  Curves 

Now  consider  the  complete  set  of  equations  of  motion  for  each  x-y-z 
component,  of  each  m-th  atom  in  the  unit  cell.  By  assuming  that  all  of  these 
motions  can  be  described  by  Hooke's  law  (and  that  each  atom  undergoes  simple 
harmonic  motion  in  response  to  any  perturbation),  the  force  on  the  m-th  atom 
in  the  unit  cell  can  be  written  as 


m 


-Mw^u 
m  m 


(52.) 


where  is  the  mass  of  the  m-th  atom,  w  is  the  frequency  of  vibration  and 
u  is  the  displacement  vector  of  the  m-th  atom  from  its  equilibrium  position. 
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as  before.  Equating  these  components  of  the  force  with  those  derived  from  the 
equations  of  motion  (Eq.  (44))  gives  a  set  of  3N  coupled  linear  equations  In 
the  terms  of  the  3N  components  of  the  atomic  displacement  vectors.  Here  N  Is 
the  number  of  atoms  In  the  unit  cell. 

There  Is  perhaps,  one  additional  complication  to  the  equations  of 

motion  that  has  been  omitted  up  to  this  point.  That  Is  the  constraint  Imposed 

by  the  periodic  boundary  conditions  associated  with  the  crystal  lattice.  These 

boundary  conditions  dictate  that  the  vector  displacements,  um  and  um, ,  of  two 

atoms  located  at  equivalent  lattice  sites  in  two  different  unit  cells  must  be 

related  by  a  phase  factor  corresponding  to  the  periodicity  of  the  lattice, 

i.e.,  u  *  =  u  exp  (ilt  •  ?  One  simple  procedure  for  assigning  the  phase 

factors  properly  is  to  multiply  each  component  of  the  displacement  vectors, 

u  ,  by  exp  (ilt  •  r  )  where  r  ,  is  the  vector  distance  from  some  fixed  origin 
m  •'  m'  m 

to  the  m-th  atom  in  the  appropriate  unit  cell. 

The  lattice  dynamical  matrix,  with  dimensions  of  3N  x  3N,  correspond¬ 
ing  to  the  right-hand  side  of  Eq.  (44)  can  be  expressed  in  terms  of  the 
coefficients  given  in  Eqs.  (46),  (50),  and  (51)  along  with  the  phase  factors 

discussed  in  the  preceedlng  paragraph.  The  matrix  element,  ,,  ,  is  defined 

11  i 

as  the  coefficient  of  the  j-th  component  of  the  k-th  displacement  vector  as¬ 
sociated  with  the  1-th  component  of  the  force  on  the  m-th  atom  in  the  unit 
cell.  It  can  be  written 


1  "i^m 
1  m 


Vkj 


TT 

m 


ijt.r_  ilt.r. 


I  ”  -  5«ke  *) 


ilt.r.  1  it.  r_ 


+  l  _L  8*mn  <C*mn1dlmi  1  '  6mke 


l  n*t 
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iit.r  iit.r 

+  D*mnfrdmnj  (5zke  “  5mke  ^ 


iit.r  iit.r 

+  E<mnidnmj'5nke  -6mke  1 


iit.r  iit.r 

F«mnid*mj(6nke  "  6mke  ) 


i  1  r_  1  r 

+  n  nL  G^m'  l  Hn>nm'idmnj  ^mke  "  6nke  ) 

n  m  *m  1  ** 


+  Imnn'idm'nj  ^6mke 


iit.r  i  it.r 

m  r  „  n 

"  6nke  J 


iit.r  i  iit.r 

+  Jmnm'idnrnj  (5m'ke  ”  {nke  ) 


iit.r,  i  it.r 

+  WiWVk*  '  finke  )1 


(53) 


where  6^  Is  the  Kronecker  delta  and  the  overall  normalization, 

exp  (-lit  •  rm)/Mn,  Is  the  mass  and  phase  factor  associated  with  Eq.  (52).  The 

various  coefficients  in  Eq.  (53)  are  defined  by: 


«m1 


*  ■(Co)*md*mi/d*m 


(54) 
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Bann  “  +^l^«mn^s1n  9imndindnm^ 


(55) 


d  ^  dnm 

Enmii  =  ^7"  dimi  cos  9imn  "  '■“imi  ‘r  anmi 
£m 


(dtmi  +  dnmi)  COS  9^  +  dnm1  COS  9^ 


(56) 


d  d 

F  .  =  d  .  +  d  -  -  d.  .  cos  9.  -  -r—  d  .  cos  6  (57) 

rimni  jmi  nmi  d  inn  imn  dnm  ntm  *mn 


G  >  =  - ( C, )„_  ^/ ( S i n2  9  -d^  d^ ,  ) 
mnm  '  1 '  mnm  '  mnm  mn  m  n 


(58) 


d2  d 

Hmnm'i  =  "T5  dmni  cos  9mnm'  *  cT^~  dm'ni  cos  6mnm' 
mn 


(59) 


d  „ 

T  «  U  *  H  *4  «  Ap  Q 

‘mnm'i  ~  *mnm'1  "  am'ni  “  dmn  umni  mnm' 


(60) 
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(61) 


29 

C4663A/bw 


g&umv+a. W»« 


Rockwell  International 

Sc  tone*  Center 


SC5266.2FR 

where  the  triple  subscript  notation  represents  one  specific  component  of  the 
vector  difference  (e.g,  d^.  =  -  ^m)j  =  (<f£)j  -  The  equations  of 
motion  can  now  be  written  in  matrix  notation  as 


N 


l 

k=l 


*  .  (u.  ) •  =  (-M  u2)(u  ). 

m.  ,k k'j  '  m  M  m'l 

*  J 


(62) 


which  has  the  form  of  an  eigenvalue  equation  where  (-M^u2)  is  the  eigenvalue  and 
(um)i  is  one  component  of  the  eigenvector.  This  eigenvalue  problem  can  be  solved 
by  numerically  diagonalizing  the  dynamical  matrix  at  any  given  wavevector  and  the 
solution  results  in  a  set  of  lattice  dispersion  curves  (u>  vs  K)  associated  with 
the  3N  eigenvalues.  In  addition,  the  3N  normal  vibrational  modes  correspond  to 
the  3N  eigenvectors  found  during  the  diagonalization  process.  The  dynamical 
matrix  elements  (Eq.  (53))  could  be  divided  by  instead  of  to  put  the 

eigenvalue  equation  in  Hermitian  form;  then  the  eigenvectors  must  be  multiplied 
by  to  obtain  the  components  of  the  actual  displacement  vectors  associated 
with  the  normal  modes.  It  should  be  emphasized  that  this  procedure  is  applicable 
to  any  arbitrary  crystal  geometry,  with  the  only  limitation  being  a  computational 
one  due  to  an  unreasonably  large  number  of  atoms  per  unit  cell. 

The  above  formalism  has  been  tested  on  GaAs,  a  zincblende  semicon¬ 
ductor.  Using  values  for  the  force  constants  of  Cg  =  47.4  eV  and  Cj  =  2.3  eV 
gives  the  dispersion  curves  shown  in  Fig.  2.  This  result  appears  to  exactly 
reproduce  Fig.  9-2  of  Reference  1.  A  more  detailed  application  was  performed  on 
the  wurtzlte  semiconductor,  CdS,  with  Cg  =  40  eV  and  Cj  =  1  eV.  First,  the 
forces  and  motion  were  assumed  to  be  constrained  to  the  x-direction  exclu¬ 
sively.  The  resulting  dispersion  curves  are  given  in  Fig.  3.  Analytical 
expressions  have  been  derived  at  K  =  0  in  order  to  check  these  curves:  us  =  0  for 
the  acoustical  mode  and 
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Fig.  2  Lattice  dispersion  curves  for  GaAs. 


31 


Rockwell  International 

Science  Center 


j 


SC5266.2FR 


^ +  fl^co +  r  ci  * 


9MjM2 


(Mj  +  M,) 


7  Cl) 


(64) 


J-  -  9C1/(M1  +  M2) 


(65) 


* 

i 

i 

j 


for  the  optical  modes.  The  appearance  of  this  low  frequency  optical  mode, 
involving  only  the  C}  force  constant,  corresponds  to  pairs  of  adjacent  planes  of 
CdS  atoms  shearing  with  respect  to  each  other.  Finally,  the  complete  dynamical 
matrix  for  CdS,  allowing  all  three  degrees  of  freedom,  was  diagonalized  for  the 
same  set  of  force  constants  (Cg  =  40  eV  and  =  1  eV)  and  the  resulting 
dispersion  curves  are  given  in  Fig.  4. 


2. 1.5. 3  Lattice  Response  to  a  dc  Electric  Field 

The  same  dynamical  matrix,  which  was  used  to  determine  the  lattice 
dispersion  curves  (Eq.  (62)),  can  be  manipulated  to  give  the  response  of  the 
crystal  lattice  to  an  external  dc  electric  field.  The  dc  field  implies  using  the 
dynamical  matrix  at  zero  wave-vector  {K  *  0).  In  order  to  find  the  perturbed 
crystal  geometry,  the  component  of  force  on  the  right-hand-side  of  Eq.  (40)  is 
replaced  with  the  i-th  component  of  the  force  exerted  by  the  dc  electric  field  on 
the  m-th  atom  in  the  unit  cell.  The  resulting  set  of  linear  equations  can  be 
written  as 


N  3 

I  I 

k*l  j*l 


e(eT>m, 


(66) 


where  again  e  is  the  electronic  charge,  (ei)  is  the  transverse  effective  charge 
corresponding  to  the  m^-th  normal  lattice  vibrational  mode  and  E^  is  the  i-th 
component  of  the  applied  dc  electric  field.  This  results  in  only  3N-3 
Independent  equations  with  3N-3  independent  unknown  rel.  •?  displacements.  By 
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Fig.  4  Lattice  dispersion  curve  for  CdS  with  general  (x,  y,  and  z) 
displacements. 
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defining  these  relative  displacements  to  be  with  respect  to  the  first  atom  in  the 
unit  cell  (k=l),  Eq.  (66)  becomes 


N  3 


l  l  \  k  (uk)  =  e(ef)m 

:=2  j=l  mi,\j  k  j  T  m1  1 


where  the  quantity  [(uk )j -(Uj )^ ]  corresponds  to  the  3N-3  unknown  relative  dis¬ 
placements.  This  above  equation  is  a  consequence  of 


®m.  k.  $m .  1 . 
k=2  i,  j  i,  j 


for  each  combination  of  m,  i  and  j  which  can  serve  as  a  convenient  computation 
check.  Equation  (67)  still  represents  a  set  of  3N  linear  equations  but  they  are 
not  all  independent;  for  example,  the  three  equations  with  m=l  can  be  expressed 
as  a  simple  linear  combination  of  the  remaining  3N-3  equations.  The  solution  to 
Eq.  (67),  which  will  yield  the  perturbed  crystal  geometry,  can  be  written 
analytically  in  determinant  notation  as 


(U1  ’j'3  *  d«CV  k.  *  (*<*?>„  E, 

'  »  J  « 


kj^.k'^j'^V  k.] 

> »  j  1 ,  j 


where  again  6k>k,  and  6jj,  are  Kronecker  delta  functions.  For  N>2,  this 
solution  becomes  unwieldy  and  a  numerical  technique  such  as  the  Gauss-Jordan 
elimination  method  Is  more  appropriate.  The  Implementation  of  this  numerical 
procedure  is  currently  in  progress  using  the  standard  and  widely  available 
mathematical  subroutine  packages.  Consequently,  results  utilizing  this 
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technique  are  not  complete  at  present  and  the  general  application  of  this 
method  Is  still  speculative. 

One  final  note  concerning  the  lattice  response  of  a  crystal  to  an 
applied  dc  electric  field  Is  a  reminder  that  the  tight-binding  formalism  for 
the  electro-optic  effect  (Eq.  (41))  depends  on  the  change  In  tk  Interatomic 
distance  with  respect  to  the  applied  dc  field.  This  change  can  be  expressed 
using  the  notation  of  Eqs.  (41)  and  (67)  as 


C(u*n),  -  (.*),] 


umni 

pdc 


(70) 


in  the  limit  of  E^c  approaching  zero  for  ^  being  the  interatomic  distance 
between  atoms  m  and  n. 


2. 1.5. 4  Transverse  Effective  Charge 


The  only  missing  ingredient  in  the  above  formalism  for  the  lattice 
response  of  a  crystal  to  an  applied  dc  electric  field  (Eqs.  (67)  and  (70)  is 
the  transverse  effective  charge,  (e*)m  .  Traditionally,  the  transverse  effec¬ 
tive  charge  has  been  introduced  in  ordlr  to  account  for  the  local  polarization 
induced  by  relative  displacements1  such  as  those  described  above  in  Eq.  (62). 
Experimentally,  there  is  a  splitting  of  the  optical  mode  frequencies  for  polar 
crystals,  even  at  zero  wave-vector  (K*0).  This  splitting14  is  a  consequence 
of  the  added  rigidity  in  longitudinal  (but  not  transverse)  lattice  modes 
resulting  from  the  local  polarization  density  induced  by  the  relative  dis¬ 
placement  of  a  pair  of  neighboring  atoms  with  unequal  charges.  Additionally, 
this  same  charge  quantifies  the  coupling  between  the  transverse  lattice 
vibrations  and  light  (i.e.,  optical  electro-magnetic  waves)  leading  to  the 
name:  transverse  effective  charge. 
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The  well-known  Lyddane-Sachs-Teller  relation  can  be  expressed  In 
terms  of  a  transverse  effective  charge  as  follows.  Let  the  optical  dielectric 
constant  be 


e»  =  1  +  4% 


(71) 


and  the  static  dielectric  constant  becomes 


e0  =  1  +  4nxe  +  4"xi  =  e»  +  4irxi 


(72) 


where  ^  is  the  electronic  susceptibility  and  xj  is  the  ionic  susceptibility. 
Now  consider  two  atoms  with  charges  of  +q  and  -q  vibrating  against  each  other, 
with  a  lattice  coupling  spring  constant,  k.  The  equation  of  motion  for  this 
system  is 


k(u2  -  Uj)  =  q(AE)  =  q(AP)/Xi 


(73) 


where  (u£  -  u^)  is  the  relative  displacement  of  the  two  atoms,  a E  is  the 
induced  change  in  the  dipole  electric  field.  The  induced  change  in 
polarization  density,  aP,  is 


&  •  l  =  q*U2  “  Ul*/n0 


(74) 
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where  Qq  is  the  volume  occupied  by  the  two  atoms.  Then,  the  ionic  suscepti¬ 
bility  becomes 


xi  =  =  P2/(no'Jl4) 


(75) 


2 

where  p  =  ^  is  the  reduced  mass.  Finally,  the  Lyddane-Sachs-Tel ler 

relation14  is 


T’  1  +  4Ve 


using  Eq.  (72).  This  can  be  rewritten  as 


(76) 


2 

*1 


2 

“r 


4nq  /  (ftnpE  ) 


(77) 


i.e.,  the  difference  in  the  square  of  the  longitudinal  and  transverse  optical 
mode  frequencies  is  proportional  to  the  effective  charge,  q,  squared. 

The  concept  of  transverse  effective  charge  can  now  be  generalized  to 
accommodate  a  crystal  of  N  atoms  per  unit  cell.  The  vibrational  motion  of 
each  atom  in  the  unit  cell  has  been  previously  described  In  Eq.  (62). 

However,  this  equation  treats  the  eigenfrequency  of  each  normal  mode  without 
accounting  for  the  change  in  local  polarization  density  Induced  by  the 
relative  displacements  of  the  atoms.  The  1-th  component  of  this  extra  force 
on  the  m-th  atom  in  the  unit  cell  is 


(78) 
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end  It  should  be  added  to  the  lattice  equation  of  motion  (l.e.,  Eq.  (62)). 
Here,  there  Is  an  effective  transverse  charge  associated  with  each  normal 
(lattice  vibrational)  mode  of  the  crystal.  Then,  if  the  approximation  that 
the  eigenvectors  for  the  longitudinal  optic  modes  are  invariant  to  this 
additional  perturbation  is  made,  a  generalized  form  of  Eq.  (77)  is  obtained: 


4TC2(e-T  )m  „ 


(79) 


where  is  the  mass  of  the  m-th  atom  in  the  unit  cell. 

For  simple  crystal  systems,  it  is  possible  to  use  Eq.  (79)  to  infer 
values  for  the  various  transverse  effective  charges  from  experimental  Lyddane- 
Sachs-Teller  splittings  in  the  optical  modes  of  the  lattice  vibrational 
spectrum.  However,  an  approach  more  aligned  with  the  philosophy  of  the  Bond 
Orbital  Model  would  be  to  develop  a  formalism  to  predict  these  transverse 
effective  charges  from  first  principles.  The  predicted  values  could  then  be 
compared  directly  to  experiment  via  Eq.  (79).  In  either  case,  the  numerical 
values  for  (e?)  could  then  be  substituted  back  into  Eq.  (67)  to  give  the 
lattice  response’of  the  crystal  to  an  external  dc  electric  field.  Finally, 
this  result  would  permit  the  numerical  evaluation  of  the  ionic  contribution  to 
the  electro-optic  susceptibility  (Eq.  (41)  to  be  completed. 

2.1.6  The  Electro-Optic  Effect  in  TeO; 

The  chemical  bonding  in  tellurium  dioxide  can  be  characterized  as  a 
mixed  tetrahedral  complex  in  which  the  tellurium  atoms  have  fourfold  coordi¬ 
nation  while  the  oxygen  atoms  have  only  twofold  coordination.  In  this  respect 
it  belongs  to  the  class  of  materials  X02,  containing  compounds  such  as  S102, 
T102,  and  AI2PO4,  which  has  received  considerable  attention.1  In  this  class 
of  complexes,  an  atom  X  is  covalently  bonded  to  four  oxygen  atoms  approxi¬ 
mately  in  a  tetrahedral  configuration.  However,  in  order  to  describe  the 
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electronic  and  lattice  features  of  crystals  from  this  class  In  a  simple  bond 
picture,  it  Is  more  appropriate  to  choose  as  a  molecular  bonding  subunit  an  X- 
0-X  structure  comprised  of  two  X  hybrid  orbitals  and  the  atomic  orbitals  of  an 
intervening  oxygen. *  Such  subunits  are  analytically  coupled  so  that  effects 
such  as  susceptibility  may  be  approximately  calculated  by  adding  the  effects 
of  independent  bonding  units.  It  is  a  generalization  of  the  Bond  Orbital 
Model  of  tetrahedral  semiconductors.  This  choice  of  bonding  subunit  greatly 
simplifies  the  calculations  in  the  relatively  complex  crystal  structure  of 
this  class  of  materials,  and  yields  a  reasonable  estimate  of  various  optical 
properties  like  the  electric  susceptibility.  Although  the  crystal  structure 
of  Te02  contains  a  large  number  of  atoms  in  the  unit  cell  and  the  detailed 
structure  of  the  Te-0  bonding  is  in  fact  more  complicated  than  that  in  Si O2 • 
similar  aspects  of  the  bonding  structure  suggest  a  parallel  approach  to  cal¬ 
culating  the  susceptibility.  Therefore  a  simplified  bond  scheme  will  be 
applied  to  Te02  in  order  to  calculate  a  specific  second-order  susceptibility: 
the  linear  dc  electro-optic  effect. 

Te02  is  an  optically  anisotropic  (positive  uniaxial)  crystal  and  is 
commonly  used  as  an  acousto-optic  material.15,16  The  linear  dc  electro-optic 
coefficient,17  other  nonlinear  optical  properties,18"22  and  elastic,  dielec¬ 
tric  and  piezoelectric  properties,23  have  been  experimentally  determined. 
Consequently,  a  comparison  of  the  Bond  Orbital  Model  prediction  for  the 
electro-optic  effect  will  be  made  with  the  experimental  result. 

2. 1.6.1  Crystal  Structure 

Tellurium  dioxide  has  the  point  group  symmetry  422,  consistent  with 
the  space  group  PA^^.24  The  lattice  constants  are  a  =  4.796A  and 
c  *  7.626A,  with  four  tellurium  atoms  and  eight  oxygen  atoms  per  unit  cell. 
Each  tellurium  has  four  oxygen  nearest  neighbors.  Each  oxygen  has  two 
tellurium  nearest  neighbors  at  different  distances  1.92A  and  2.09A  while  the 
Te-O-Te  bond  angle  is  140.8°.  The  properties  considered  here  vary  as  the 
square  of  the  distance  and  an  average  of  2.00  A  will  be  used.  The  crystal¬ 
lographic  unit  cell  is  shown  schematically  in  Fig.  5. 
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Fig.  5  TeC^  unit  cell  with  nearest  neighbor  bonds:  (a)  30  perspective, 
(b)  viewed  down  the  x-axIs,  and  (c)  viewed  down  the  z-axls. 
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2. 1.6. 2  First-Order  Susceptibility 

Following  the  Bond  Orbital  Model  formalism  presented  for  Si O2 • ^  the 
linear  susceptibility  of  Te02  can  be  obtained  by  calculating  the  polariz¬ 
ability  of  each  Te-O-Te  molecular  unit  and  then  summing  the  eight  resulting 
dipole  moments  over  the  unit  cell  (normalized  by  the  volume  of  the  unit  cell). 

As  an  approximation,  the  electronic  structure  of  a  symmetric  Te-O-Te 
subunit  may  be  determined  by  considering  an  atomic  p-state  on  each  tellurium 
atom  to  be  spatially  directed  toward  the  oxygen  atom.  This  atomic  p-state 
must,  in  the  real  system,  be  some  sort  of  s-p  hybridized  orbital  since  each 
tellurium  atom  has  four  nearest  neighbor  oxygen  atoms  while  only  three 
p-states  are  available  for  these  four  bonds.  The  two  unoccupied  tellurium  p- 
states  in  the  Te-O-Te  subunit  will  only  significantly  interact  with  the  two 
occupied  oxygen  p-states  in  the  x-z  plane,  as  Illustrated  in  Fig.  6a.  (The 
oxygen  s-state  is  ignored  because  its  energy  is  considerably  lower  than  any 
tellurium  state.  The  s-state  on  the  tellurium  Is  not  Included  since  It  is 
occupied  and  the  only  unoccupied  states  are  orthogonal  tellurium  p-states. 
Lastly,  the  py-state  on  the  oxygen  is  also  orthogonal  to  these  tellurium 
s-states).  To  further  simplify  the  calculation,  only  the  ppa  coupling  will  be 
included  since  it  is  a  factor  of  four  larger  than  the  ppn  contribution. 
Referring  to  the  states  Identified  in  Fig.  2a,  the  bonding  (Bx  and  Bz)  and 
anti  bonding  (Ax  and  Az)  states  for  the  symmetric  Te-O-Te  unit  in  the  x  and  z 
directions,  respectively,  are: 
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T«  #1 


Te  #2 


Fig.  6  Te02  molecule  subunit:  (a)  atomic  orbital  assignments  and 
(b)  coordinate  system  for  unsymmetric  subunit. 
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Bx  =  ?  ^  "  8x  +  |P2>)  +  ^  +  2BX  | Px >3 

Bz  =  \  [•'l  -  3Z  (- |Pj>  +  |p2>)  +  &  +  2SX  |PZ>] 

(80) 

Ax  =  7  +  ~ex  ( |px>  +  |p2>)  -  ^  -  2ex  |px>] 

Az  ■  7  C^r  +  sz  (-|Pi>  -  |p2>)  +  J2  -  2BZ  |pz>] 

where 


»*  =  Vj//*$  ♦  ^ppa  s,"‘  ♦ 

(81) 

‘  V/V3  *  2VpZpa  ‘°s‘  ♦  • 

The  polar  energy  for  Te02  is 
ETe  -  e° 

V3  =  P-  2~  P  =  3.59  eV  (82) 

while  the  Vppo  plays  the  role  of  a  covalent  energy  for  the  symmetric  subunit. 
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ppo 


-2.22 


mdf 


=  4.22  eV 


(83) 


using  the  parameters  introduced  in  Ref.  (11). 

To  calculate  the  polarizabilities  of  the  Te-O-Te  subunit,  the 
perturbation  of  the  polar  energy,  in  response  to  the  optical  electric  field, 
is  considered.  This  perturbation  is  eExx  or  eEzz  for  the  x  and  z  directions, 
respectively.  Writing  out  the  x  and  z  bonding  wavefunctions  in  a  first-order 
perturbation-theory  expansion  of  the  antibonding  wavefunction  gives  a 
polarizability  of 


2e2  <Bi|xi|Ai>2 
%  • 


(84) 


where  and  e ^  are  the  bonding  state  and  anti  bonding  state  energies  given 
by 


(85) 


Evaluating  the  dipole  matrix  elements  In  the  usual  fashion,1  the 
polarizabilities  become 


°x  =  —  dLS31,n2t  tBZ(1  -  BZ2)  +  M1  •  BX2)]  (86) 


2  e2  d2  cos2*  ”  0*0z^  0x0z 
'  — * —  c~ - 
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Numerical  evaluation  of  Eq.  (86)  for  Tet^  yields 


ax  =  +0.533  A3 
az  =  +5.04  A3 


(87) 


The  first-order  susceptibility  tensor  component. 


is  defined  by 


P  s  £ 

’  }  ^ 


£j 


(88) 


where  is  the  induced  polarization  density  and  Ej  is  the  optical  electric 
field.  We  could  obtain  a  susceptibility  by  averaging  the  effect  of  each  unit 
over  angle,  giving  (c^  +  az)/3  (since  Oy  =  0  in  our  description)  and  multiply 
by  the  density  of  bonding  units  N0  =  8/(a^c)  (or  oxygen  atoms): 


,(D 


=  N„ 


az)/3 


0.085. 


(89) 


As  another  alternative,  the  tensor  Xi j  for  the  real  crystal  is  determined  from 
the  bonding  subunit  polarizabilities  by  summing  over  the  eight  subunits  in  the 
unit  cell  as  follows.  The  electric  field  along  the  j-th  crystal  coordinate 
axis  must  be  transformed  Into  an  electric  field  vector  In  the  coordinate 
system  of  each  Te-O-Te  subunit.  The  x  and  z  components  of  this  transformed 
electric  field  induce  a  dipole  moment 


P 


Exx  + 


Ez2 


(90) 
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in  the  Te-O-Te  subunit  coordinate  system.  This  dipole  moment  then  must  be 
transformed  back  into  the  crystal  coordinate  system.  Finally,  xjj^  is 
calculated  by  summing  up  the  i-th  components  of  the  transformed  dipole  moments 
from  each  of  the  eight  Te-O-Te  subunits  and  dividing  by  the  volume  of  the  unit 
cell.  Note  that  the  transformation  between  the  crystal  axes  and  each  Te-O-Te 
subunit  coordinate  system  is  different  (see  Table  1). 

In  a  coordinate  system  of  the  crystallographic  axes,  the  calculated 
linear  susceptibilities  are 


Y(i>.  YU)  - 
X11  x22 


0.092 


0.071 


(91) 


Note  that  these  first-order  susceptibilities  are  simply  the  sum  of  the 

Individual  bond  polarizabilities  and  do  not  include  any  corrections  due  to 

bond-bond  interactions  or  local  field  effects.  Corrections  of  this  type  are 

sometimes  semiempi rically  taken  into  account  by  introdyuing  a  scaling  factor, 

Y1  This  scaling  factor  enters  in  as  y^  and  it  is  used  to  match  the 

f  1 1  1 J 

calculated  values  of  x^  with  those  derived  from  the  experimentally  measured 
refractive  indices  using  the  relation 


'V  i + 


(92) 


Comparison  with  the  experimental  values  for  the  ordinary  and  extraordinary 
Indices,  nQ  =  2.26  and  ne  *  2.41  respectively,18  yields  an  average  y  =  2.1. 
The  significance  of  this  value  Is  that  it  can  now  be  used  for  the  calculation 
of  the  second-order  susceptibility  without  introducing  additional  scaling 
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Table  1 

The  Transformation  Between  the  Crystal  Axes  and  the  Te-O-Te  Subunit 
Is  Specified  by  the  Direction  Cosine  of  Each  Crystal  Axis 
With  Respect  to  Each  Subunit  Axis.  The  Direction  Cosines  Are: 

Crystal  Coordinates  Crystal  Coordinates  Crystal  Coordinates  Crystal  Coordinates 
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parameters.  It  should  be  pointed  out  that  the  calculated  sign  of  the  bire¬ 
fringence,  an  =  ng  -  nQ,  is  predicted  Incorrectly.  Within  the  context  of  the 
Bond  Orbital  Model,  the  sign  of  the  birefringence  is  specified  strictly  by 
crystal  geometry.  These  differences  in  the  ordinary  and  extraordinary  indices 
are  perhaps  be  due  to  the  bond-bond  coupling  or  anisotropic  local  field 
effects. 

2. 1.6.3  Second-Order  Susceptibility 

The  second-order  susceptibility,  x^j»  may  be  divided  into  two 
parts:  an  electronic  part  and  an  ionic  displacive  one.  The  electronic 
contributions  arise  from  the  changes  In  the  first-order  electronic  polariz¬ 
ability  due  to  the  application  of  a  dc  electric  field  while  keeping  all  atomic 
positions  fixed.  Calculations  in  Te02  (point  group  symmetry  422)  show  that  this 
electronic  contribution  vanishes  in  accordance  with  the  Kleinman's  symmetry 
criterion.2  The  ionic  displacive  contribution  is  determined  by  perturbing  the 
crystal  geometry  with  a  dc  electric  field,  as  described  below,  and  then  recal¬ 
culating  the  first-order  susceptibility  using  the  formalism  given  in  Sec¬ 
tion  2. 1.6. 2.  We  hold  the  crystal  surfaces  rigid  or  clamped  and  then  do  not 
Include  any  distortion  arising  from  the  piezoelectric  effect.  The  second- 
order  susceptibility,  and  hence  the  "clamped"4  electro-optic  coefficient,  is 
then  proportional  to  the  difference  between  the  two  linear  susceptibilities 
associated  with  the  perturbed  and  unperturbed  geometries. 

In  determining  the  ionic  displacive  contribution  in  the  clamped  case, 
the  only  important  degrees  of  freedom  for  the  electro-optic  effect  correspond 
to  the  IR-active  optical  modes  in  which  the  displacement  of  an  atom  gives  rise 
to  a  net  change  In  the  electrical  dipole  moment.  Since  the  ratio  of  the  Te 
mass  to  the  0  mass  Is  about  eight,  one  can.  In  principle,  decouple  the  problem 
of  finding  the  acoustical  modes  from  that  of  the  optical  modes.  Therefore, 
only  the  vibrations  of  the  oxygen  atoms  in  an  immobile  tellurium  lattice  will 
be  considered  and  any  coupling  between  the  various  local  oxygen  modes  will  be 
neglected.  (Note:  a  detailed  description  of  the  vibrational  modes  in 


48 

C4663A/bw 


Rockwell  International 

Science  Center 


SC5266.2FR 

paratellurite  has  been  given  elsewhere.^)  Thus,  again  consider  a  bonding 
subunit  Te-O-Te  as  in  Fig.  6  and  now  distinguish  the  two  bond  lengths  which 
were  neglected  In  the  preceding.  Letting  the  elastic  motion  of  the  oxygen  be 
characterized  by  a  "bond-stretching"  force  constant,  C0,  which  describes  the 
change  In  Te-0  interatomic  distance,  and  a  "bond-rocking"  force  constant,  Cj, 
which  describes  the  tilt  of  the  Te-O-Te  unit  out  of  the  x-z  plane,  the  elastic 
energy  can  be  written 


7  Co 


(d, 


V 


1 

7 


(d; 


'  d2}  +  1  c 

12  +  7  C1 


(6S)‘ 


(93) 


where  d'j  and  d 1 ^  are  the  time- dependent  Te-0  bonding  distances,  d^  and  d2  are 
the  corresponding  equilibrium  distances  and  <$£  is  the  angular  tilt.  Here,  the 
"bond- bending"  force  for  the  Te-O-Te  bond  angle  has  been  neglected  because  it 
Is  much  smaller  than  the  "bond-stretching"  force  and  any  oxygen  displacement 
in  the  x-z  plane  causes  changes  in  bond  length.  F or  the  coordinate  system 
shown  in  Fig.  6b,  the  equations  of  motion  for  the  oxygen  atom  (to  first-order 
in  the  displacement)  are 


-  M  u,2a*  =  -  C0[(-J-  +  -J-)  Ax  -  (■ 


X  Z.  X  z, 

'?  4 


C1 

-  Maj2  Ay  =  -  v  Ay 
xo 


(94) 


-  M(u2Az  = 


Cnt-  tr  +  -T2-)  4x  +  (4-  +  %-)  42] 


J1 


where  Zj  and  Z2  specify  the  tellurium  positions  and  x0  gives  the  location  of 
the  oxygen  atom.  (Note  that  the  motion  In  the  y-dlrection  is  completely 
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uncoupled  from  the  x  and  z  motions.)  In  the  approximation  of  a  symmetric 
Te-O-Te  unit  where  dj  =  62  and  zi  =  z2»  three  normal  vibrational  modes  are 
purely  along  the  x,  y  and  z  axes  with  the  eigenfrequencles  given  by 


2  C  sin 29 

a2 - 2 - 

*  M  d2 


C1 

2  _  1 
y "  m  x2 

0 


“z2  = 


2  c  cos2e 


M  d2 


(95) 


respectively,  where  0  =  (*  -  e')/2  and  9*  is  the  Te-O-Te  bond  angle.  For  the 
unsymmetrlc  case,  the  equations  of  motion  can  also  be  obtained  in  closed  form 
and  the  resultant  set  of  vibrational  modes  and  frequencies  only  deviate  a 
small  amount  from  the  symmetric  case. 

The  three  local  mode  frequencies  u»y  and  1^  will  be  the  same  for 
all  the  oxygen  atoms  in  the  unit  cell.  Allowing  for  some  small  frequency 
distribution  for  the  various  oxygen  atoms  due  to  coupling  between  different 
modes  through  the  tellurium  atoms  in  the  real  crystal,  this  local  mode 
approximation  leads  to  three  sets  (or  bands)  of  optical  frequencies,  each 
centered  at  the  local  mode  frequency  given  by  Eq.  (95).  Indeed  the  experi¬ 
mental  Infrared  reflectivity  data26  show  distinct  bands  which  can  be  inter¬ 
preted  as  arising  from  the  oxygen  local  modes.  This  predicts  Immediately  a 
ratio  =  tan  0  =  0.36.  Assigning  the  centers  of  the  two  highest  bands  to 

the  frequencies  of  u>®xp  t  =  0.70  x  1014  sec-1  and  *  1.45  x  1014  sec”1 

for  the  bond-stretching  local  modes  and  the  third  one  with  a  frequency 
of  u>®xp  *  *  0.33  x  1014  sec”1  for  the  rocking  mode  of  the  Te-O-Te  subunit,  we 
find  =  0.48.  The  following  estimates  for  C0  and  Cj  are  obtained  from 

Eq.  (95): 
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C0  =  110  eV 
Cj  *  8.0  eV 


(96) 


where  C0  was  determined  from  the  average  fit  to  u>xxp  *  and  w®xp  *. 

The  magnitude  of  these  values  for  C0  and  Cj  are  in  reasonable  agreement  with 
estimates  for  bond-stretching  and  bond-bending  force  constants  in  semicon¬ 
ductors.*  It  is  worth  mentioning  here  that  due  to  the  decoupling  of  the  bond¬ 
stretching  modes  and  above)  from  the  Te-O-Te  bond-rocking  mode  (u»y), 
the  frequency  of  the  latter  Is  determined  by  only  and  is  therefore  "softer" 
than  the  bond-stretching  modes.  As  it  turns  out,  this  softer  mode  makes  the 
major  contribution  to  the  electro-optic  coefficient  since  it  is  most  affected 
by  the  dc  electric  field. 

In  order  to  estimate  the  motion  of  the  oxygen  atoms  in  response  to  an 
applied  dc  electric  field,  we  introduce  a  transverse  effective  oxygen  charge 
eT*  associated  with  each  mode.1  It  Is  defined  to  be  the  electric  dipole 
Interaction  by  a  single  oxygen,  divided  by  the  di splacement.  Expressions  to 
determine  these  effective  charges  in  terms  of  the  Bond  Orbital  Model 
parameters  already  Introduced  have  previously  been  derived1  for  the  Si -0-Si 
unit.  Using  the  same  expressions  for  Te-O-Te: 


(ef)y  =  Z*  =  Bx  +  Bz  (97) 

(ef)x  =  Z*  +  (3  sin2e  -1)  Bx(l  -  B2)  +  3  s1n2e  Bz(l  -  B2) 

(e*)z  -  l*  +  (3  cos2e  -  1)  Sz(l  -B2)  +  3  cos29  Bx(l  -  Bx) 
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where  0X  and  6Z  are  given  in  Eq.  (81).  Numerically  evaluating  the  parameters 
for  a  symmetric  Te-O-Te  subunit  gives 


(e*)  =  1.40 

'  T  '  X 


(ef)  =  1.42 


(98) 


(e*)z  =  2.60. 


Experimentally,  these  tranverse  effective  charges  give  the  local 
polarization  induced  by  the  relative  oxygen  displacements  and  hence  the 
splitting  between  the  transverse  and  longitudinal  optical  modes.  For  each 
mode,  the  transverse  effective  charge  is  shown,  in  Section  2. 1.6. 5  below,  to 


be  given  by  a  Lyddane-Sachs-Tel ler14  relation 


“L0  "  “TO 


4u(e*) 2  e2/3mOoEj 


(99) 


where  e  is  the  electronic  charge,  M  is  the  mass  of  the  oxygen,  is  the 
average  volume  per  oxygen  atom  and  e^(=  ej  is  the  optical  dielectric  con¬ 
stant.  Using  this  relation  along  with  the  approximate  experimental  splittings 
2  2  2 

of  (ui^q  -  Wj-q)  =  .093,  .028  and  .326  sec"41  for  the  x,  y  and  z  modes  respec¬ 
tively,  we  obtain  experimental  transverse  effective  charges  of 


Rockwell  International 

Sc  tone*  Center 


SC5266.2FR 


(ef)x  -  1.75 

(®j)y  -  0.96  (100) 

(e*)z  =  3.27 


which  can  be  directly  compared  with  the  theoretical  estimates  given  in 
Eq.  (98). 

The  force,  ,  in  the  i-th  direction  on  an  oxygen  atom  is  given  by 


F 


i 


e(e*). 


(101) 


where  E?c  is  i-th  component  of  the  dc  electric 
of  the  Te-O-Te  subunit.  Finally,  substituting 
the  components  of  the  displacements,  which  can 


field  in  the  coordinate  system 
this  force  into  Eq.  (94)  gives 
simply  be  written  as 


AX 


i 


M  w* 


(102) 


for  the  symmetric  Te-O-Te  subunit.  (The  displacements  for  the  unsymmetric 
subunit  can  also  be  obtained  In  an  exact  (but  more  complicated)  closed  form.) 
These  displacements  can  now  be  used  to  generate  a  perturbed  crystal  geometry 
for  an  appropriate  dc  electric  field  and  the  first-order  susceptibility  can  be 
recalculated,  as  In  Section  2. 1.6. 2. 
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The  defining  relation  for  the  second-order  (electro-optic) 
susceptibility,  x^]»  Is 


P.  =  S(XW 


l  xjj?  Ef )  E°pts  l  L 
k  1Kj  K  J  j 


E0pt 


(103) 


where,  as  before,  P^  is  the  induced  polarization  density,  xij^  is  the  first- 

order  susceptibility,  E?pt  is  the  optical  electric  field,  e!/  is  the  dc  elec- 
T  J  x 

trie  field  and  x^  represents  the  "total"  susceptibility.  The  electro-optic 

susceptibility  tensor  is  then  obtained  from  the  difference: 


4*L  ■  xlj(EdC'°)  -  xlj(EaC-0)  ■  l  Eac 


(104) 


where  each  particular  component  of  the  tensor  is  determined  by  considering  a 
specific  direction  for  ^c. 


The  procedure  for  predicting  this  electro-optic  susceptibility  is 
somewhat  similar  to  that  used  for  the  unperturbed  linear  susceptibility. 
First,  a  dc  electric  field,  E^c,  is  applied  in  the  crystal  coordinate  system 
of  each  Te-O-Te  subunit  and  the  relative  vector  displacement  of  each  oxygen, 
induced  by  this  dc  field,  is  calculated  according  to  Eq.  (102).  Next,  the 
vector  displacement  of  each  oxygen,  given  in  the  Te-0-Te  subunit  coordinate 
system,  must  be  transformed  back  into  the  crystalline  axes  thus  determining 
the  perturbed  crystal  geometry.  The  new  linear  susceptibility  can  now  be 
calculated  and  the  electro-optic  susceptibility  obtained  from  the  difference. 
The  process  was  carried  out  numerically  for  fields  in  each  of  the  three 
crystal  directions  for  E?c.  One  computational  difficulty  was  encountered 

*  j.  j. 

involving  the  magnitude  of  e£  .  If  Ej;  was  too  small,  then  the  change  in 
linear  susceptibility  was  less  than  the  computer  round-off  error  in  the 
calculation.  On  the  other  hand,  if  e£c  was  too  large,  third-order  effects 
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began  to  appear.  Realizing  this  problem  and  choosing  an  appropriate  value  for 
the  magnitude  for  E^c  gave  a  3  *  3  x  3  tensor  for  the  electro-optic 
susceptibility,  in  which  all  components  are  zero  except  for 


x213  =  x312  =  “x123  =  ”x321  =  0,37  *  10  ?  statvolt  ^105^ 


consistent  with  the  symmetry  of  point  group  422.  In  this  numerical 
evaluation,  the  result  listed  in  Eq.  (105)  has  been  scaled  by  y2  as  discussed 
in  Section  2. 1.6. 2.  Finally,  the  clamped  electro-optic  coefficient  is  given 
by 


-  4tt 


r41 


TT 

e  0 


^13 


=  -0.15  x  10 


-7 


cm 

statvolt 


-  -  0.51  x 


10”*2m/v 


(106) 


where  the  values  for  ne  and  nQ  are  theoretical  values  calculated  from  Eq.  (92) 
with  the  linear  susceptibility  also  scaled  using  the  same  y. 

An  alternative  to  using  the  Te-0-Te  subunits  for  calculating  the 
linear  and  nonlinear  susceptibilities  Is  to  employ  the  bond-by-bond  approach, 
as  has  been  typically  done  with  the  tetrahedral  semiconductors.**  There,  the 
first-order  susceptibility  can  be  written 


e2 

TV 


l 

bonds 


(V 


T 


(107) 


where  e  is  the  electronic  charge,  y,  V2  and  V3  are  defined  in  terms  of  sp^ 
hybridization  In  Ref.  11,  dm  Is  the  m-th  component  of  the  Interatomic  distance 
vector  between  the  bonding  atoms,  v  is  the  volume  of  the  unit  cell  and  the 
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summation  Is  over  all  bonds  in  the  unit  cell.  Applying  this  method  to  the 
sixteen  bonds  in  the  Te02  unit  cell  gives  first-order  susceptibilities  of 


(108) 


which  can  be  compared  with  Eq.  (91)  after  accounting  for  the  use  of  a  slightly 
different  y  -  1.77  as  listed  in  Ref.  12. 

The  electro-optic  susceptibility  and  electro-optic  coefficient  can 
then  be  determined  within  the  context  of  this  bond-by-bond  method  by  using 
Eqs.  (103),  (104)  and  (105).  The  results  are 


(109) 


(110) 


which  compare  reasonably  well  with  Eqs.  (105)  and  (106)  even  though  the 
assumptions  about  tetrahedral  coordinate  that  went  Into  Eq.  (107)  are  somewhat 
dubious  for  Te02. 

2. 1.6.4  Discussion 

The  relatively  close  agreement  In  the  results  of  the  molecular 
Te-0-Te  subunit  approach  and  the  bond-by-bond  approach  suggests  that  the 
details  of  the  electronic  structure  may  not  be  as  Important  as  the  crystal 
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geometry.  This  is  because  the  non-zero  tensor  elements  of  the  susceptibility 
arise  when  the  components  of  the  dipole  moments  do  not  exactly  cancel  when 
summed  over  the  unit  cell.  The  extent  of  the  cancellation  Is  strictly  a 
geometrical  effect  and  the  results  are  relatively  Insensitive  to  the 
parameters  that  represent  the  electronic  structure. 

An  experimental  value  for  the  "unclamped"  or  constant  stress  electro¬ 
optic  coefficient,  r^.  In  Te0£  has  previously*7  been  determined  to  be 


rjj  =  -  0.76  x  10‘12  m/V 


(111) 


The  above  calculations  predicted  the  "clamped"  or  constant  strain 
coefficient,  r^,  which  can  be  related  to  the  "undamped"  coefficient  by5 

r41  =  r41  “  p44  d14  (112) 


where  P44  is  a  shear  photoelastic  constant  and  d^  Is  the  piezoelectric 
constant.  The  magnitudes  of  both  P44  and  d^  have  previously  been 
measured:*5*27  |p^|  =  0.17  and  1  |  =  0.81  x  10”11  C/N  . 

Substituting  In  these  numbers  Into  Eq.  (112)  yields  two  possible 
experimental  values  for  the  "unclamped"  electro-optic  coefficient,  depending 
on  the  sign  of  the  product  (p44  d^).  For  (p44  dj4)  negative. 


rjj  =  +  0.62  x  10"12  m/V 
and  for  (P44  d14)  positive. 


(113) 
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rjj  =  -  2.14  x  10"12  m/V  . 


(114) 


This  pair  of  numbers  can  now  be  directly  compared  with  Cqs.  (106)  or  (110). 

If  (P44  d^)  is  less  than  zero,  then  the  magnitude  of  the  predicted  electro¬ 
optic  coefficient  agrees  remarkably  well  with  experiment  but  a  discrepancy 
exists  in  the  sign.  The  theory  determines  this  sign  solely  from  the  crystal 
geometry.  On  the  other  hand,  if  (p44  d^)  is  greater  than  zero,  then  the  sign 
agrees  but  the  magnitude  of  the  prediction  is  approximately  a  factor  of  four 
too  small.  The  theory  predicts,  then,  that  p44  and  d14  have  the  same  sign. 

2. 1.6. 5  The  Lyddane  Sachs  Teller  Relation  For  TeO? 

The  optical  vibrational  modes  for  Te02  can  be  rather  well  separated 
into  three  distinct  sets,  with  frequencies  uy,  uy,  and  a^.  If  there  were  no 
dipoles  associated  with  them,  these  oscillators  would  be  treated  as  identical, 
with  N0  modes  per  unit  volume,  where  Nq  is  the  number  of  oxygen  atoms  per  unit 
volume. 

To  understand  the  small  shifts  in  frequencies  due  to  electrostatic 
forces,  consider  a  particular  set  of  oscillators,  say  those  of  frequency, 
uy.  These  are  degenerate  vibrational  modes,  one  for  each  oxygen,  and  could  be 
equivalently  written  as  modes  with  well-deflned^v^umber^y  giving  each 
local  oscillator  a  displacement  amplitude  of  ue  J  2  ,  where  r^  is  the 
atomic  position,  and  this  will  In  fact  give  the  modes  which  interact  with 
infrared  optical  fields  of  a  corresponding  wavenumber. 

As  a  first  problem,  assume  that  every  oscillator  moves  in  a  direction 
exactly  parallel  to  the  wave-vector,  £,  which  has  been  chosen.  (This  will 
never  be  the  case,  but  will  make  the  calculation  clear  for  the  interesting 
case.)  If  the  dipole  moment  Is  given  by  (e^)z  times  the  displacement  of  the 
atom  (this  Is  the  definition  of  the  transverse  charge  (e|)z),  then  the 
polarization  density  is  P  -  N0(e|)z  u  e^k,rj  ”  where  u  is  the  vector 
amplitude  parallel  to  £.  The  electric  field  Is  obtained  from  the  charge 
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accumulation,  -i(c*  through  Poisson’s  Equation.  By  dividing  by  the 

electronic  dielectric  constant,  the  net  field,  and  hence  the  force  on  each  atom 
given  by  e(e|)z  times  the  electric  field.  Is  determined  and  leads  to  a  shift 
in  frequency.  This  gives  the  famllar  Lyddane-Sachs-Teller  result  for  the 
difference  between  the  longitudinal  frequency,  ojj_z ,  and  the  transverse 
frequency,  u>yz,  (for  which  there  Is  no  charge  accumulation  and  therefore  no 
shift  in  frequency).  It  is 

4  -  <4 =  4i,<e?>z  e2vMc-  •  (n5> 


If  we  now  examine  the  relevant  case  for  TeOg  where  the  orientation  of 

the  polarization  vector  for  the  mode  is  different  for  each  different 

oxygen,  similar  charge-polarization-modulated  modes  result  which,  without 

coupling,  all  have  the  same  frequency.  However,  In  this  case,  the  dipole  of 

each  mode  makes  an  angle  0j  with  ic  so  Its  contribution  to  the  polarization 

density  Is  reduced  by  a  factor  of  cos  0j.  Furthermore,  the  force  felt  by  the 

dipole  along  its  direction  of  motion  due  to  the  field  of  all  other  dipoles  Is 

reduced  by  cos  0j  so  It  behaves  exactly  as  if  It  were  aligned  with  the  field 

but  had  Its  effective  charge  reduced  by  a  factor  of  cos  0j.  If  every  dipole 

had  the  same  0^,  the  Lyddane-Sachs-Teller  relation  given  above  would  be  valid 
oJ  9  9 

but  with  (e|)z  replaced  by  (ef)z  cos^Oj.  Clearly  In  our  averaged  treatment  of 
the  local  modes,  cos20j  should  be  re-laced  by  Its  average,  one  third,  for  the 
case  of  random  orientations. 

The  distinction  between  longitudinal  and  transverse  modes  in  this 
scheme  must  be  clarified.  The  phase  for  each  local  mode  was  chosen  so  that 
the  sign  of  the  polarization  along  £  was  the  same  for  neighboring  atoms, 
corresponding  to  a  longitudinal  mode.  Had  a  particular  transverse  direction 
been  selected  and  the  sign  In  that  direction  chosen  to  be  the  same  for  neigh¬ 
boring  atoms,  a  mode  with  transverse  polarization  would  have  been  obtained  but 
the  longitudinal  components  would  have  cancelled  on  average,  and  there  would 
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have  been  no  shift  In  frequency.  Clearly,  at  long  wavelengths,  there  are 
longitudinal  and  transverse  optical  modes  just  as  there  are  longitudinal  and 
transverse  acoustic  modes.  There  is  no  difficulty  for  the  long  wavelengths  of 
Interest  here  and  the  relation  for  our  system  becomes 


■L  -  4  *  <ef>z  eV<3Me->  • 


(116) 


The  corresponding  expressions,  with  the  same  factor  of  one  third,  may  be 
written  for  the  x-  and  y-modes. 

2.2  Experimental 


2.2.1  The  Electro-Optic  Effect  in  Tl-^  As  Se-* 

We  have  completed  measurements  of  the  linear  electro-optic  tensor 
component  (rc  =  +1.5  x  lO^m/V)  for  the  chalcogenide  salt  TljAsSej  (TAS).  TAS 
crystallizes28  the  trigonal  space  group  R3m,  with  the  point  group  3m.  This 
implies  the  electro-optic  tensor  has  the  following  form:^ 


r  n 

0 

r13 

-rll 

0 

r13 

f  0 

0 

r33 

r1j = 

0 

r51 

0 

1 

r51 

0 

0 

0 

-rll 

0  / 

(117) 


The  dc  electric  field  was  applied  along  the  crystal lographnlc  ^  axis.  With 


60 

C4663A/bw 


Rockwell  International 

Science  Center 


SC5266.2FR 

the  light  propagating  along  polarized  45°  from  £  towards  and  the  dc 
field  applied  along  _z,  the  dielectric  Impermeability  tensor*  has  the  form 


—  +  r  E 
2  r13t3 


2  +  r33E3 


Multiplying  by  the  appropriate  polarization  tensors,  we  find: 


pBp=l-(0  11)  0  4-+rnE. 

ft  nf  13  3 


_ _ +  r  E 

2  33t3 


‘J(V  r13E3  +  T2  +  r33E3 ) 


=  B„  +  AB 
0 
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which  Implies  that 


r33  ^ E  3 


(120) 


For  this  propagation  direction,  the  effective  index  is: 


ne  = 


2nn 
e  o 


/2"(n^  +  njj 


(121) 


If  we  assume  the  incident  wave  is  linearly  polarized  at  45°  to  the 
optical  axis  upon  entry  into  the  crystal,  the  propagating  wave  can  be 
decomposed  into  the  two  orthogonal  modes  (ordinary  and  extraordinary)  in  the 
crystal.  As  the  wave  travels  through  the  crystal,  these  two  modes  get  out  of 
phase  with  one  another  due  to  their  different  propagation  velocities.  The 
accumulated  phase  difference  r  at  the  exit  surface  of  the  crystal  is  given  by: 


r  = 


no> 


(122) 


where  i  =  the  crystal  length  and  \Q  =  the  free-space  wavelength.  Thus  the 
emerging  beam  will,  in  general,  be  elliptically  polarized  due  to  the  arbitrary 
phase  between  the  energy  In  each  of  the  two  modes.  Experimentally,  we  may 
quantify  the  elllpticity  through  the  use  of  a  Babinet  Compensator  and  a  final 
linear  polarizer  analyser.  The  application  of  a  transverse  dc  electric  field 
may  affect  this  propagating  beam  In  three  ways.  First  because  TAS  Is  also 
piezoelectric,  the  electric  field  may  change  the  physical  dimensions  of  the 
sample  due  to  electrically  Induced  mechanical  strain: 
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In  this  case  =  E3  so  we  have 
el  =  d31E3 

e2  =  d31E3  (124) 

e3  =  d33E3 

The  two  strain  components  ancl  e3  d0  not  produce  a  direct  mechanical  effect 
on  the  sample  which  effects  our  measurement.  However  ej  which  Is  a  contrac¬ 
tion  (or  expansion)  of  the  crystal  along  _x,  will  change  the  physical  length  of 
the  sample.  This  in  turn  will  produce  a  change  in  the  phase  at  the  exit 
surface.  Although  no  direct  experimental  measurement  of  the  piezoelectric 
effect  in  TAS  exists,  we  can  assume  that  it  is  small  consistent  with  similar 
semiconducting  materials  like  GaAs.  Thus  the  direct  mechanical  perturbation 
will  be  ignored  and  the  measurement  will  provide  only  the  "unclamped"  value 
for  the  electro-optic  coefficient.  There  Is,  however,  a  second-order  indirect 
effect  often  referred  to  as  the  converse  piezo-optic  effect  In  which  the 
electrically  Induced  strain  couples  to  the  refractive  index  via  the  strain- 
optic  effect.  The  strain-optic  effect  produces  on  index  change: 


An 


1  *  '  5“  p1jej 


(125) 


where  p^j  Is  the  strain-optic  coefficient  and  cj  is  the  strain  component  as 
before.  Then 


An. 


-  f  pudkjEk 


(126) 
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or  we  may  define  an  effective  converse-piezoelectric  electro-optic  tensor 
component : 


1jdkj 


(127) 


The  photo-elastic  tensor 
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Again  taking  E 
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and 
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elements  of  both 
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(128) 


and  the  d^  tensor,  we  find  the  remaining  non-zero  components  are: 

4 

{An)l  "  ‘  T  ^ plld31  +  p12d32  +  p13d33^E3 
_3 


2  ^pll  +  p12^d31  +  p13d33^E3 


(129) 
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(An)2  =  -  j-  (P2id3i  +  P22d32  +  p23d33^E3 


“  T~  ^p12  +  p22^d31  +  p13d33^E3 


(An)3  =  -  j-  (P31d31  +  P32d32  +  p33d33^E3 


-  j~  ^2^p31d31 ^  +  p33d33^E3 


(130) 


Thus  we  note  that  both  the  and  components  of  the  Index  (ordinary  and 
extraordinary  respectively)  will  be  affected  by  the  dc  elecrlc  field.  At 
present,  data  for  the  magnitude  of  pjj  and  p33  have  been  taken  at  this 
laboratory  under  company  sponsored  IRSD  program,  but  no  data  exists  for  the 
remaining  photoelastic  components  or  the  piezoelectric  tensor  components. 

Thus  although,  as  shown  below,  we  were  able  to  measure  a  value  for  the 
"electro-optic"  effect.  It  is  at  present  not  possible  to  quantitatively 
separate  out  effect  due  to  the  piezoelectric  effect  and  the  converse  piezo¬ 
optic  effect. 

2.2.2  Measurement  of  the  Electro-Oplc  Effect  In  TAS 

A  single  crystal  boule  of  TAS  was  grown  from  the  melt  using  the 
vertical  Bridgman-Stockbarger  technique.29  An  experimental  sample  of  dimen¬ 
sions  1.27  cm  x  0.54  cm  x  0.56  cm  was  wire  sawed  from  the  boule,  aligned  using 
back  reflection  Laue  photographs,  and  mechanically  polished.  The  sample  geom¬ 
etry  Is  as  shown  In  Fig.  7.  In  our  early  experiments  a  wire  grid  polarizer 
was  used  as  an  analyzer  after  the  crystal.  A  compensator  was  placed  between 
the  sample  and  the  analyzer.  The  compensator  and  analyzer  were  then  adjusted 
so  that  essentially  no  light  was  transmitted  through  the  analyzer  to  the 
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Fig.  7  Tl2AsSe3  crystal  orientation  for  rc  electro-optic  measurement 
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detector.  As  the  dc  electric  field  was  slowly  Increased,  a  series  of  oscil¬ 
lations  In  the  signal  level  at  the  detector  ^as  observed  (see  Fig.  8).  How¬ 
ever,  further  experiments  showed  that  these  preliminary  results  were  the 
result  of  joule  heating  of  the  sample  due  to  an  excessive  current  flowing 
through  the  sample  when  high  voltage  was  applied.  Indeed  It  was  possible  to 
reproduce  the  oscillating  transmission  by  blowing  hot  air  from  an  electric 
heat  gun  across  the  sample  with  no  applied  dc  voltage.  This  excessive  current 
was  produced  by  a  low  value  for  the  bulk  resistivity  in  the  sample.  Rather 
than  pursue  the  traditional  methods  to  remove  Impurities  and  defects  to  obtain 
Intrinsic  resistance,  which  tend  to  be  extremely  costly  and  time  consuming,  we 
have  chosen  an  altrnatlve  approach.  A  high  voltage  pulser,27  shown  In  Fig.  9, 
was  used  to  apply  pulses  of  high  voltage  to  the  sample.  The  pulse  width  was 
kept  short  and  the  repetition  rate  correspondingly  low  so  that  the  duty  cycle 
was  reduced  and  consequently  heating  effects  were  reduced  if  not  eliminated. 
The  average  power  dissipated  In  the  sample  was  typically  of  the  order  of 
1  mW.  A  high-speed  (150  MHz)  i-urrent-to-voltage  amplifier  was  used  In 
conjunction  with  a  small  area  InSb  photodiode  to  observe  the  change  In  optical 
transmission  produced  by  the  applied  high  voltage  pulse.  It  should  be  pointed 
out,  however,  as  shown  In  Fig.  10,  the  pulse  was  sufficiently  long  as  to  allow 
time  for  mechanical  response  of  the  sample.  That  Is  to  say,  the  measurement 
was  made  under  constant  stress  conditions  (unclamped).  The  signal  from  the 
amplifier  was  then  fed  to  one  channel  of  an  oscilloscope  and  displayed  simul¬ 
taneously  with  the  high  voltage  pulse.  The  process  of  taking  data  consisted 
of  (1)  setting  the  amplitude  of  the  high  voltage  pulse,  (2)  reading  the  ampli¬ 
tude  of  the  optical  signal  as  a  function  of  analyzer  angle  (over  90®  of  rota¬ 
tion)  and  (3)  then  repeating  this  process  at  several  peak  pulse  voltages.  One 
can  show  that  for  propagation  of  light  polarized  at  an  angle,  a,  relative  to 
the  optical  axis,  through  a  blrefrlngent  medium  and  through  an  analyzing 
polarizer  oriented  at  an  angle,  e,  (relative  to  the  optical  axis),  the 
resulting  transmission  will  be: 
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Intensity  modulation  as  a  function  of  applied  dc  electric  field  In 
Tl3AsSe3. 
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Fig.  10  High  voltage  pulses  applied  to  TloAsSe-j  and  resulting  modulation  of 
analyzed  light  Intensity  transmitted  through  Tl3AsSe3. 
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If  the  propagating  medium  is  an  electro-optic  material  oriented  similar  to  our 
sample,  then  an  applied  dc  electric  field  will  modify  the  term  containing  an 
via  the  electro-optic  effect.  The  difference  in  Intensities  with  the  electric 
field  on  and  off  is 


al(0,a,E) 


sina  cosa  sine  cose 


(to 


"0 
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eff 


(132) 


This  change  in  transmitted  Intensity  as  a  function  of  analyzer  angle,  0,  for  a 

given  "dc"  electric  field  Is  precisely  the  data  we  have  obtained 

experimentally.  The  data  can  then  be  numerically  fitted  to  the  above 

eff 

expression  to  yield  a  value  for  rc  .  Figure  11  shows  a  fit  to  the 
experimental  data.  This  process  yields  a  value  of  r®ff  «  +  1.5  x  1012  m/V  at 
T  =  300K  and  \  *  3.39  un.  As  previously  Indicated,  more  data  Is  needed  to 
explicitly  solve  for  each  of  the  various  contributions  to  this  value.  Future 
studies  on  TAS  hopefully  will  supply  this  information. 
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3.0  CONCLUSION 

In  summary,  the  microscopic  understanding  of  the  electro-optic  effect 
has  progressed  substantially  during  the  course  of  this  program.  A  tight- 
binding  theory  for  dielectric  susceptibilities  has  been  specifically  applied 
to  the  linear  dc  (Pockel's)  effect.  The  distinction  between  electronic  and 
ionic  contributions  was  explicitly  maintained.  A  formalism  to  address  the 
generalized  lattice  dynamic  problem,  associated  with  the  Ionic  contribution  to 
the  electro-optic  effect,  has  been  developed  using  a  method  which  should  be 
applicable  to  all  materials,  e.g.,  ferroelectric  crystals  that  exhibit  large 
electro-optic  coefficients. 
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